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ABSTRACT 


Laser-Doppler  velocimeter  measurements  were  made  in  the  flow  field  of  a  stalled 
cascade  of  controlled-diflRision  stator  blades.  Tests  were  conducted  at  10  degrees  of 
incidence  above  the  design  inlet-flow  angle  in  order  to  verify  previous  measurements.  A 
unique  method  of  data  presentation  was  offered  to  characterize  the  unsteady  positive  and 
negative  velocities  about  their  mean.  Laser-sheet  flow  visualization  was  performed  and 
showed  the  unsteadiness  of  the  stalling  phenomena  within  the  blade  row  passages. 
Additionafly,  a  leading-edge  separation  bubble  was  observed  and  reversed  flow  was 
measured  within  the  bubble  using  the  laser-Doppler  anemometer.  Power-spectral  density 
and  auto-correlation  analysis  results  are  presented.  In  parallel  with  the  experimental 
measurements,  a  computational  fluid-dynamics  study  was  initiated  in  an  attempt  ultimately 
to  predict  stall.  Viscous  flow  computations  were  completed  at  on-design  and  near-stall 
conditions.  Pressure  distributions,  separation  bubble  re-attachment  locations,  and  loss 
predictions  were  compared  to  previous  experimental  studies.  Generally  good  agreement 
was  obtained  for  on-design  conditions  with  degradation  in  prediction  nearing  stall. 
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I.  INTRODUCTION 


Compressor  stall  and  oflf-design  behavior  impose  performance  constraints  on  gas 
turbine  engines.  The  open  literature  currently  lacks  adequate  documentation  of 
quantitative  measurements  of  turbomachinery  flow  fields  at  and  beyond  stall.  This 
knowledge  is  essential  for  predicting  the  onset  of  stall  with  computational  methods  and 
optimizing  blade  geometries  for  off-design  range  with  minimal  degradation  in  design-point 
performance.  Towards  optimum  design-point  performance,  controUed-dififtision  (CD) 
blading  was  developed.  The  purpose  of  these  blades  was  to  avoid  boundary  layer 
separation  and,  in  the  case  of  supercritical  (SC)  blading,  to  be  shock-fi-ee  in  the  transonic 
range. 

Extensive  laser-Doppler  velocimeter  (LDV)  investigations  [Refs.  1  -  3]  of  flow 
through  a  controlled-diflEusion  compressor  cascade  have  been  conducted  in  the  region 
below  stall  at  the  Naval  Postgraduate  School  (NFS).  Recently,  the  reported  cascade  was 
stalled  with  an  inlet-flow  angle  of  50  degrees  and  an  extensive  mapping  of  the  flow  field 
was  performed  using  a  two-component  LDV  system  by  Ganaim  [Ref  4].  The  motivation 
for  the  present  study  was  to  partially  map  the  flow  field  near  stall  and  to  establish  the 
repeatability  of  the  results  reported  in  Reference  4.  Between  Ganaim’s  work  and  the 
present,  the  tunnel  was  reconfigured  tAvice  to  pre-stall  flow  angles.  This  was  done  to 
confirm  the  repeatability  of  the  flow  field  and  to  expand  the  available  database  for 
imsteady  separated  flows.  Also,  since  the  previous  data  were  presented  on  a  mean-flow 
basis,  it  was  desired  to  more  properly  characterize  the  mean  flow  into  its  positive  and 
negative  components.  This  required  a  special  method  of  post-processing  the  velocity- 
histogram  data,  as  described  within  this  report. 

It  was  also  desired  to  characterize  the  unsteady  flow  using  power-spectrum  and 
auto-correlation  analyses.  It  was  hypothesized  that  there  could  exist  a  fundamental 
fi-equency  for  the  stalled  cascade,  that  it  could  be  determined  fi-om  the  measurements  made 
in  the  present  study  and  should  match  that  of  the  previous  data  [Ref  4]. 

In  parallel  with  the  LDV  study,  a  computational  fluid  dynamics  (CFD)  study  was 
initiated  in  an  attempt  ultimately  to  predict  stall.  Previous  CFD  studies  on  the  present 
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cascade  include  those  of  Elazar  [Ref  1],  Stow,  et  al  [Ref  5],  Hobson  [Ref  6],  and  Kang, 
et  al  [Ref  7]. 

Before  attempting  to  calculate  the  unsteady,  stalled,  flow  through  a  cascade  (as 
was  measured  in  the  present  study),  it  was  necessary  first  to  evaluate  the  code’s  predictive 
ability  at  incidence  angles  ofF-design  and  near  the  stall  point.  A  quasi-three-dimensional 
computer  code  [Refs.  8-9]  for  the  solution  of  the  thin-layer  Navier-Stokes  equations  was 
used  to  predict  the  cascade  performance  from  design  (40  degrees  inlet-flow  angle)  to  near 
stall  (48  degrees  inlet-flow  angle).  The  code  was  used  to  successfully  compute  the  blade 
surface  pressure  distribution,  the  leading  edge  separation  bubble  and  total  pressure  loss 
coefficient,  for  inlet  flow  angles  fi'om  40  to  48  degrees. 
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11.  EXPERIMENTAL  SETUP 


A.  LOW-SPEED  CASCADE  WIND  TUNNEL 

The  subsonic  cascade  wind  tunnel  contained  20  controlled-diffiision  blades.  Each 
blade  was  10  inches  in  span,  5.01  inches  in  chord  and  separated  by  3  inches  in  the 
pitchwise  direction.  The  cascade  solidity  (blade  chord-to-pitch  ratio)  was  1.67  and 
Reynolds  number  based  on  chord  was  approximately  700,000.  A  detailed  description  of 
the  test  facility,  test  section  and  CD  blading  was  fully  documented  in  Reference  1.  Figure 
1  shows  a  schematic  of  the  cascade.  The  tunnel  was  configured  for  an  inlet-flow  angle 
(3i)  of  50  degrees,  and  an  exit-flow  angle  (^2)  such  that  the  downstream  end-wall  static 
pressure  distribution  was  uniform. 

B.  INSTRUMENTATION 

Blade  surface  pressure  measurements  were  recorded  using  a  48-channel  Scanivalve 
system  controlled  by  an  HP-9000  computer.  Figure  1  shows  the  “fully-instrumented 
blade”  location  within  the  test  section.  Reference  1  gives  the  location  of  the  blade  surface 
pressure  points.  A  full  description  of  the  data  acquisition  system  is  given  in  Reference  10. 

A  two-component,  four-beam  TSI  Model  9100-7  LDV  system  was  used  for  all 
flow  field  measurements.  The  LDV  setup,  including  laser  type,  optics,  atomizer,  seeding, 
traverse  mechanism  and  data  acquisition  were  thoroughly  documented  in  Reference  4. 

The  flow  was  visualized  by  injecting  fog  particles  of  2  pm  in  size  through  a  wand 
inserted  through  the  end  wall  upstream  of  the  blading.  A  300  mW  air-cooled  argon-ion 
laser  was  used  as  the  light  source  for  laser-sheet  flow  visualization.  The  laser  beam  was 
transmitted  to  a  lens  system  via  a  fiber-optic  cable.  The  lens  system  ended  with  a 
cylindrical-piano  lens,  which  produced  the  laser  sheet.  The  laser  sheet  was  introduced  into 
the  cascade  tunnel,  ahead  of  the  blades,  through  the  sidewall.  A  schematic  of  the  system, 
showing  the  location  of  the  introduction  of  the  fog  and  the  position  of  the  video  camera 
(located  outside  the  window),  are  shown  in  Figure  2. 
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LOCATION 
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Instrumented 


Figure  2.  Laser-Sheet  Flow-Visualization  System 
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c. 


EXPERIMENTAL  PROCEDURE 


1.  Surface  Pressure  Measurement 

The  tunnel  plenum  chamber  pressure  was  set  at  12  inches  of  water  (gauge)  and  the 
blade  surface  pressure  measurements  were  taken  as  described  in  Appendix  B  of  Reference 
10.  Inlet  total  and  static  pressures  were  recorded  2  1/3  chord  lengths  ahead  of  the  blades. 
These  pressures  were  used  to  non-dimensionalize  the  blade-surface  pressure 
measurements  as  a  coefficient  of  pressure;  i.e.,  Cp  =  (piocai  -  Poo)  /  (ptoo  -  Poo)- 

2.  LDV  Measurements 

A  full  listing  of  the  experimental  runs  is  contained  in  Appendix  A.  Blade-to-blade 
surveys  were  conducted  at  eleven  stations  ahead  of,  between,  and  downstream  of  the 
blade  row.  All  surveys  were  taken  by  collecting  3,000  data  points  at  each  location  in  the 
survey.  All  LDV  data  were  collected  in  the  “coincidence  mode”  with  a  2.0x10*  psec 
window  used  for  validation.  Table  1  summarizes  the  survey  stations  and  locations.  Figure 
3  illustrates  the  same  information  graphically.  The  numbers  in  parentheses  on  the  figure 
indicate  the  horizontal  distance  between  measurement  points.  The  present  study  gave  an 
average  inlet-flow  angle  of  50.05  degrees.  Reference  4  results  were  for  an  average  Pi  of 
50.21  degrees. 


Station 

Region 

Location  (%  Axiid  Chord) 

1 

Inlet 

30.9  %  forward  of  LE  (1.5  inches  forward  of  LE) 

lb 

(C 

4.28  %  “ 

Id 

1.00  %  “ 

2a 

Passage  (Suction  Side) 

0%  @LE 

3 

u 

5.15%  aftofLE 

7 

a 

30.9  %  “ 

(1.5  inches  aft  of  LE)  I 

14 

u 

93.6  %  “ 

15 

M 

98.7  %  “ 

17 

Wake 

6.18%  aflofTE 

18 

cc 

12.7  %  “ 

19 

j 

20.6  %  “ 

(1.0  inches  aft  of  TE)  1 

Table  1.  LDV  Survey  Stations  and  Locations 
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2 


Figure  3.  LDV  Survey 


(inches) 


;ions  Between  Blades  6  and  7 
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Additionally,  eleven  specific  locations  were  measured,  with  an  increased  number  of 
data  points,  for  power-spectrum  and  auto-correlation  analyses.  Table  2  below  summarizes 
the  measurements  and  the  number  of  data  points  gathered  for  each  point.  Figure  4 
depicts  the  stations  at  which  these  dedicated  runs  were  performed.  The  points  close  to  the 
blade  were  chosen  to  be  within  the  intermittent  reverse-flow  region,  and  the  outboard 
points  were  chosen  to  be  outside  the  reverse-flow  region. 


'  Locaticm 


(0.0000,  -6.2920) 


(-1.1331,-4.7905) 


(-0.8394, 4.5420) 


(-0.1225,  -3.2920) 


(0.3988,  -1.7921) 


(1.0000,  -1.7921) 


(0.2500, 0.3620) 


(1.5000, 0.3620) 


(0.2500,  1.0620) 


(1.2500,  1.0620) 


Number  of  Points  Acquired 


3K,  30K,  50K 


3K,  30K _ 


3K,  30K _ 


3K,  30K 


3K,30K 


30K 


30K 


30K 


30K 


30K 


30K 


Table  2.  LDV  Power-Spectrum  /  Auto-Correlation  Locations 


The  tunnel  reference  velocity  (Vref)  for  each  run  was  computed  as  described  by 
Elazar  in  Reference  1.  The  results  are  shown  in  Appendix  B.  The  reference  velocity  was 
used  to  non-dimensionalize  the  data  for  comparison  with  previous  studies. 

Prior  to  each  run,  the  probe  volume  was  always  referenced  to  the  same  location 
between  the  blades  using  an  alignment  tool.  This  allowed  for  a  high  confidence  level  in 
specifying  station  position  for  comparison  with  previous  studies.  The  alignment  procedure 
was  developed  and  described  in  detail  in  Reference  1 . 
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Figure  4.  LDV  Power-Spectrum  /  Auto-Correlation  Locations 
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a.  Inlet  Surveys 

The  inlet  region  was  surveyed  across  stations  1,1b,  and  Id  (Figure  3).  The 
survey  spanned  three  and  one-third  passage  widths  in  the  pitchwise  direction  (-9  inches  to 
+1  inch)  for  station  1,  and  one  and  one-third  passage  widths  in  the  pitchwise  direction  (-2 
inches  to  +  2  inches)  for  stations  lb  and  Id.  The  laser  beam  system  was  horizontal  when 
measuring  at  station  1.  A  4  degree  pitch-up  of  the  LDV  table  was  required  for  stations 
lb  and  Id  to  avoid  blade  shadowing  of  the  vertical  beam,  and  to  achieve  closer  access  to 
the  leading  edge.  The  potential  problems  due  to  pitching  and  yawing  the  LDV  system 
were  treated  by  Hobson  and  Shreeve  in  Reference  3.  Their  analysis  showed  the  maximum 
spatial  error  due  to  probe  volume  orientation  to  be  0.3  mm,  which  is  the  minimum 
diameter  of  the  measurement  volume. 

The  optics  were  configured  the  same  throughout  the  present  study  with  the 
488-nm  blue  beam  measuring  the  horizontal  (U)  component  and  the  514.5-nm  green  beam 
measuring  the  vertical  (V)  component.  Frequency  shifting  was  performed  as  outlined  in 
Reference  4  in  order  to  capture  the  large  reversed-flow  velocities. 

b.  Passage  Surveys 

The  passage  surveys  were  performed  between  blades  6  and  7.  All  passage 
surveys  extended  1.8  inches  (approximately  60  %  passage  width)  away  from  the  suction- 
side  surface  of  blade  6.  The  LDV  system  was  pitched  up  2  degrees  and  yawed  4  degrees 
left  for  stations  2a,  3  and  7  for  better  access  to  the  suction  side  of  the  blade.  Stations  10, 
14  and  ITwere  measured  with  only  4  degrees  of  left  yaw. 

c.  Wake  Surv^s 

The  LDV  system  was  reset  to  zero  pitch  and  zero  yaw  for  measurements  at 
stations  17,  18,  and  19.  The  wake  surveys  were  centered  around  blades  6  and  7  and 
spaimed  two  passage  widths  (6  inches). 

d.  Laser  Power-Spectrum /Auto-Correlation  Measurements 

It  was  desired  to  expand  the  analysis  of  several  points  throughout  the  flow 
field  by  collecting  a  larger  amount  of  data  points.  Instead  of  the  standard  3,000  points 
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gathered  at  each  location,  30,000  data  points  were  chosen  for  the  locations  listed  in  Table 
2.  This  choice  corresponded  to  the  limiting  number  of  data  points  that  could  be  acquired 
at  positions  with  low  data  rates  such  that  the  acquisition  software  reached  its  maximum 
timeout  of  30  minutes.  This  corresponded  to  a  minimum  data  rate  of  16.7  Hz  in  order  to 
collect  30,000  samples.  Stations  2a  and  3,  near  the  separation-bubble  region,  were  most 
affected  by  this  limitation. 

3.  Flow  Visualization 

The  laser-sheet  flow  visualization  was  photographed  through  the  optical  access 
window  using  a  hand-held  8mm  camera  at  30  fi’ames/second  and  a  VHS  camera  at  60 
fi’ames/second.  The  filming  was  performed  at  night  to  reduce  glare  and  optimize  the  flow 
field  image  contrast.  A  total  of  approximately  15  minutes  of  footage  was  obtained  of  the 
intermittent  reversed  flow  inside  the  blade  passage  and  of  the  flow  near  to  the  passage 
entrance,  including  the  leading-edge  separation  bubble. 
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III.  COMPUTATIONAL  FLUID  DYNAMICS 


A.  GRID  GENERATION 

The  grid  generation  program  employed  was  a  modified  version  of  GRAPE  (GRids 
about  Aerofoils  using  Poisson’s  Equation),  [Ref  11],  This  version  allowed  for  the 
generation  of  a  periodic  C-type  grid  for  turbomachinery  cascades.  The  namelist  input  file, 
along  with  a  description  of  the  effort  put  forth  in  refining  the  leading  edge  geometry  for 
the  grid,  is  described  in  Appendix  E. 

The  grid,  a  portion  of  which  is  depicted  in  Figure  5,  contained  240  grid  points 
wrapped  around  the  blade  in  the  streamwise  direction  and  49  grid  points  in  the  blade-to- 
blade  direction.  The  grid  was  much  finer  in  the  normal  direction  near  the  surface  for 
improved  viscous  flow  resolution,  particularly  of  the  leading  edge  separation  bubble.  The 
outer,  or  periodic,  boundaries  corresponded  in  size  to  the  cascade  wind  tunnel  as 
described  in  Experimental  Setup.  The  inlet  region  was  located  30%  of  an  axial  chord 
upstream  of  the  leading  edge.  The  wake  region  extended  two  chord  lengths  downstream 
of  the  trailing  edge. 
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Figure  5.  Two-Dimensional  Periodic  C-Grid 
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B.  COMPUTATIONAL  FLOW  SOLVER 


Numerical  computations  were  made  using  Rotor  Viscous  Code  Quasi-3D 
(RVCQ3D)  of  Reference  8.  RVCQ3D  solved  the  thin-layer  Navier-Stokes  equations  by 
employing  an  explicit,  multi-stage,  Runge-Kutta  algorithm.  It  could  account  for  the 
effects  of  rotation,  radius  change  and  stream  surface  thickness  variation  in  blade-to-blade 
turbomachinery  flows.  The  Baldwin-Lomax  algebraic  eddy- viscosity  model  [Ref  12]  was 
used  to  characterize  the  turbulence  within  the  flow  field. 

Computations  were  performed  at  the  following  inlet  flow  angles:  40,  43,  46,  and 
48  degrees.  This  was  done  to  allow  comparisons  with  existing  experimental  results. 
Validation  of  the  computed  results  was  achieved  by  comparing  pressure  coefficient 
distributions  around  the  blade  surface,  leading-edge  separation  bubble  sizes,  and  loss 
computations. 

A  typical  namelist  input  file  is  shown  in  Appendix  F.  All  solutions  were  run  to 
6,000  iterations  and  axial-velocity-density  ratio  (AVDR)  was  held  at  6%  for  all 
computations.  The  AVDR  was  modeled  by  defining  two  stream  surface  points  at  the  inlet 
and  exit  using  the  namelist  6  parameter  ‘NMN’  as  shown  in  Appendix  F. 

Although  the  inlet  flow  angle  was  specified  initially,  the  flow  solver  computed  to  a 
new  inlet  condition  based  on  the  prescribed  exit-static  to  inlet-stagnation  pressure  ratio. 
The  inlet  flow  angle  varied  due  to  changes,  primarily,  in  the  tangential  velocity 
component,  as  the  axial  velocity  did  not  change  due  to  continuity  constraints. 

Solutions  were  computed  using  the  National  Aerodynamic  Simulation  (NAS) 
Program  CRAY  C90  supercomputer  and  the  NPS  CRAY  Y-MP  EL98  supercomputer. 
CPU  run  time  for  the  NAS  C90  was  between  3-4  minutes  and  the  NPS  EL98  required 
about  55  minutes.  Analysis  of  the  solution  files  was  carried  out  using  FAST  graphics 
software  [Ref  13]  and  separate  FORTRAN  programs  given  in  Appendices  H  (blade- 
surface  pressure  distribution)  and  I  (computed  inlet-region  flow  angle). 
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IV.  RESULTS  AND  DISCUSSION 


A,  BLADE-SURFACE  PRESSURE  COEFFICIENT 

A  comparison  of  pressure  coefficient  around  the  blade  surface  between  the  present 
study  and  Reference  4  is  shown  in  Figure  6.  The  results,  after  reconfiguring  the  tunnel, 
showed  reasonable  agreement  with  the  pressure  distributions  obtained  earlier  at  the 
leading  edge;  particularly  on  the  pressure  side.  There  was  some  difference  near  the 
trailing  edge  which  was  possibly  due  to  a  difference  in  setting  of  the  tunnel’s  exit-plane 
tailboards.  The  setting  of  the  tailboards  significantly  affected  the  exit  static-pressure 
distribution. 


Figure  6.  Pressure  Distribution  Comparison 
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B.  LDV  SURVEY  COMPARISONS 

The  data  analysis  and  comparisons  which  follow  were  performed  without  any 
refinement  bounds  for  the  velocity  histograms.  All  the  results  from  Reference  4  were  re¬ 
analyzed  in  this  manner. 

1.  Inlet  Survey  (Station  1) 

Figures  7  and  8  show  the  non-dimensionalized  mean  velocity  and  turbulent 
intensities,  respectively.  The  two  studies  showed  variations  in  the  mean  velocity  of  less 
than  2%  across  three  passage  widths.  The  period  in  the  variation  in  the  mean  velocity 
corresponded  to  the  3  inch  blade  separation,  indicating  that  the  blades  were  felt  30%  of  an 
axial  chord  ahead  of  the  leading  edge.  The  U  and  V  turbulence-intensity  levels  were 
similar,  and  this  indicated  the  turbulence  was  isotropic  at  the  inlet.  The  increase  in 
turbulence  intensity  from  2%  to  4%  over  the  extent  of  the  survey  is  noted.  This  behavior 
was  not  observed  at  lower  inlet-flow  angles. 

2.  Passage  Survey  (Station  7) 

Mean-flow  results  for  station  7  are  presented  in  Figure  9.  There  was  good 
agreement  across  the  passage  which  further  established  repeatability  with  Reference  4’s 
results.  The  turbulence  intensity  in  Figure  10  also  agreed  well  and  showed  anisotropic 
behavior  of  the  turbulence  near  the  blade  surface  out  to  50  %  pitch.  Beyond  this  point  the 
turbulence  became  isotropic  in  the  inviscid  region,  as  the  levels  of  turbulence  in  the  U  and 
V  components  were  comparable.  The  data  were  presented  as  turbulence  intensities,  yet 
because  the  flow  field  was  intermittently  reversible  the  standard  deviation  of  the 
histograms,  given  by  the  LDV,  did  not  represent  the  turbulence  in  the  flow.  Rather,  these 
should  be  looked  upon  as  a  measure  of  the  unsteadiness  in  the  flow,  which  was  a 
combination  of  the  random  turbulence  and  intermittent  reverse  flow. 
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Figure  7.  Velocity  Comparison  at  Station  1 


Station  1  Turbulence  Intensity 


Figure  8.  Turbulence  Intensity  Comparison  at  Station  1 
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Figure  9.  Velocity  Comparison  at  Station  7 
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Figure  10.  Turbulence  Intensity  Comparison  at  Station  7 
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Figures  1 1  and  12  separate  the  mean  horizontal  (U)  and  mean  vertical  (V)  velocity 
components  into  their  respective  positive  and  negative  mean  values.  The  velocity 
histograms  for  each  point  across  the  passage  were  specially  edited  such  that  only  negative 
velocities  were  retained  and  averaged  to  create  a  negative  velocity  trace  and,  likewise, 
only  positive  velocities  were  retained  and  averaged  at  each  point  to  create  a  positive 
velocity  trace.  It  was  apparent  from  Figures  1 1  and  12  that  the  weighted  average  of  the 
negative  trace  and  positive  trace  would  make  up  the  mean  velocity  for  each  component 
across  the  passage  as  previously  shown  in  Figure  9. 

This  analysis  showed  more  of  the  character  of  the  intermittent  reversed  flow  than 
did  a  plot  of  the  mean  velocity.  What  occurred  in  the  unsteady  stalling  process  was  not 
represented  by  the  arithmetic  mean.  The  envelope  about  the  mean  contained  the  positive 
and  negative  velocities  that  truly  characterized  the  flow’s  behavior.  Since  the  data  were 
taken  in  coincidence  mode,  one  disadvantage  of  this  form  of  presentation  was  that  during 
editing  of  the  positive  U  velocity-component  data,  all  of  its  negative  points  were 
discarded,  as  well  as  the  corresponding  V  velocity-component  data.  The  opposite  but 
equivalent  process  occurred  during  editing  of  the  V  velocity  component;  such  that  all  the 
data  that  were  left  were  the  velocity  components  in  the  first  quadrant  of  a  flow  that  was 
oscillating  with  components  in  all  four  quadrants.  It  was  apparent  from  the  figures  that  by 
30  %  pitch,  the  flow  no  longer  exhibited  reverse  velocities  in  either  the  U  or  V 
components. 
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Figure  11.  U  Velocity  Comparison  at  Station  7 
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Figure  12.  V  Velocity  Comparison  at  Station  7 
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3.  Wake  Survey  (Station  19) 

The  wake  surveys  were  compared  at  station  19.  Figures  13  and  14  show  the 
results  and  indicate  good  general  agreement  for  the  mean  velocity  and  turbulence 
intensity/unsteadiness  across  two  passage  widths  for  the  two  studies.  The  periodic 
influence  of  the  blades  was  strongly  evident  in  the  V  component  of  Figure  13.  Also,  the 
asymmetry  of  the  wake  was  especially  apparent  with  the  larger  velocity  deficit  on  the 
suction  side  of  the  blade.  Figure  14  once  again  points  out  the  anisotropic  behavior  of 
turbulence/unsteadiness  in  this  same  region. 

Figures  15  and  16  are  the  results  of  the  special  editing  as  outlined  for  station  7  in 
Figures  1 1  and  12.  It  was  of  interest  to  note  that  the  V-component  reversed-flow 
velocities  at  station  19  were  of  the  same  magnitude  as  the  V-component  reversed-flow 
velocities  at  station  7. 
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Figure  13.  Velocity  Comparison  at  Station  19 
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Figure  14.  Turbulence  Intensity  Comparison  at  Station  19 
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Figure  15.  U  Velocity  Comparison  at  Station  19 


Figure  16.  V  Velocity  Comparison  at  Station  19 


C.  PASSAGE  SURVEY  (STATION  3) 

Measurements  obtained  at  32  locations  at  station  3  are  shown  in  Figures  17-20.  A 
detailed  analysis  to  resolve  the  unsteady  characteristics  of  the  flow  was  made  at  the  three 
locations  (1,  7  and  20)  shown  in  Figure  17  (U  component  of  velocity)  and  Figure  19  (V 
component  of  velocity).  The  histograms  of  the  U  and  V  components  of  velocity  at  the 
three  locations  are  shown  in  Figure  18  and  Figure  20,  respectively.  Point  1  was  located 
just  off  the  blade  surface.  This  was  inside  the  leading-edge  separation  bubble,  evidenced 
by  the  negative  mean  velocity  of  all  the  data  points  in  the  histogram  (Figures  18a  and  20a). 
The  second  location  (point  7)  was  at  2.5  %  pitch.  This  point  was  within  the  intermittent 
reversed-flow  region,  as  indicated  by  the  bi-modal  velocity  histograms  (Figures  18b  and 
20b).  The  third  location  (point  20)  was  at  17  %  pitch,  and  was  outside  the  region  of 
intermittent  reversed-flow. 

Figure  17  shows  the  mean  pitchwise  velocity  along  with  the  separate  averages  of 
the  positive  and  negative  components.  The  progression  from  a  negative  mean  velocity  to 
a  positive  mean  velocity  was  apparent  with  the  crossover  distinctly  evident  in  the  bi-modal 
histogram  at  point  7.  It  was  of  interest  to  note  that  the  arithmetic  mean  for  the  flow  rarely 
occurred.  The  negative  and  positive  velocity  peaks  were  predominant,  which  was 
indicative  of  the  unsteady  stalled  flow. 

The  streamwise  and  pitchwise  velocity  component  behavior  at  station  3  appear  to 
be  similar  because  the  blade  surface  was  oriented  at  approximately  45  degrees  to  the  two 
pairs  of  laser  beams  (X-Y  axes  for  the  velocity  measurements).  Figure  21  provides  an 
alternate  view  of  the  previous  data  by  plotting  the  pitchwise  (U)  and  axial  (V)  velocity 
components  together. 
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Figure  17.  U  Velocity  Survey  at  Station  3 
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Figure  18.  U  Velocity  Histograms  at  Station  3 
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Figure  19.  V  Velocity  Survey  at  Station  3 
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Figure  20.  V  Velocity  Histograms  at  Station  3 


D.  POWER-SPECTRUM  /  AUTO-CORRELATION  ANALYSIS 

Power  spectral  density  and  auto-correlation  analyses  were  performed  on  the  data 
collected  as  previously  indicated  in  Table  2.  Only  station  3  (point  7)  and  station  7  (points 
3  and  20)  have  been  included  in  Appendk  C.  Therein  are  contained  the  velocity 
histograms  and  time-history  plots  for  these  points.  This  analysis  was  performed  using 
TSI’s  FIND  software,  version  4.03.  The  settings  chosen  for  the  analysis  were  the 
following;  Window  Type:  Hcaming^  Fast  Fourier  Transform  (FFT)  Method:  Random 
Correlation,  and  Correlation  Type:  Random  Correlation.  In  all  cases  the  default  lagged- 
time  slot  (Ax)  was  used  in  order  to  satisfy  the  Nyquist  sampling  criterion. 

It  was  initially  believed  that  a  large  data  set  was  required  to  properly  carry  out  this 
type  of  analysis.  The  data  sets  acquired  for  station  3,  point  7,  included  3K  and  3  OK  data 
points.  The  comparison  of  the  two  sets  showed  similarity  in  their  behavior.  The  same  can 
be  said  for  station  7,  points  3  and  20,  which  were  compared  between  3K,  3  OK,  and 
Ganaim’s  3K  data  points.  The  only  point  of  recurring  significance  throughout  the  power 
spectrum  and  auto-correlation  plots  was  in  the  region  between  20-25  Hz  as  shown  in 
Figure  22  (30K  points,  present  study),  and  Figure  23  (3K  points.  Reference  4)  for  station 
7,  point  3.  No  further  interpretation  of  the  results  obtained  was  attempted  due  to  the  many 
limitations  of  discretely-sampled,  randomly-arriving  data  points.  Additionally,  since  this 
routine  involved  averaging  the  data  over  the  lagged-time  slot  (Ax),  there  was  some 
question  whether  the  true  character  of  the  flow  had  been  captured. 
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Figure  22.  Power  Spectrum  for  30K  Points 
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Graph  Component  J 
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Figure  23.  Power  Spectrum  for  3K  Points 


The  use  of  larger  data  sets  did  not  appear  to  be  of  significance  in  the  final  analysis. 
The  primary  limitation  turned  out  to  be  the  data  sampling  rate.  The  maximum  frequency 
for  spectral  analysis  was  limited  to  1/2  the  data  rate  due  to  the  Nyquist  criterion,  which 
prevents  aliasing  of  the  data.  Table  3  shows  the  data  rates  and  maximum  fi-equencies  for 
the  spectral  analysis  plots  of  Appendix  C. 


Location 

Number  of  Data  Pointe 

Data  Rate 
(Hz) 

Maximum 

Frequency 

(Hz) 

Stations,  Point? 

SK 

S9 

19 

a 

SOK 

49 

24 

Station?,  Points 

SK 

82 

41 

(6 

SOK 

12? 

6S 

cc 

SK  (Ganaim) 

606 

SOS 

Station?,  Point 20 

SK 

69 

34 

(C 

SOK 

?1 

S5 

SK  (Ganaim) 

61S 

S06 

Table  3.  Spectral-Analysis  Maximum  Frequencies 


E.  FLOW  VISUALIZATION 

Visualization  of  reversed-flow  velocities  within  a  blade  passage  was  documented 
on  VHS  videotape.  The  tunnel  was  run  at  a  plenum  chamber  pressure  of  12  inches  H2O  to 
match  the  conditions  during  the  LDV  measurements.  The  shutter  speeds  for  the  Smm 
recording,  were  varied  fi’om  l/250th  to  1/1 000th  and  l/4000th  of  a  second  with  the  rate  of 
30  fi'ames/second  remaining  constant.  It  was  found  that  the  slower  shutter  speed  of 
l/250th  of  a  second  was  better  able  to  capture  the  streaklines  of  the  intermittent  reversed- 
flow  due  to  longer  particle  traces  for  each  image. 

Upon  replaying  the  videos  at  normal  and  slower  speeds,  several  observations  were 
made  for  both  the  present  study  and  the  flow  visualization  performed  in  Reference  4.  The 
reversed-flow  region  appeared  to  be  confined  to  the  blade  passage  and  did  not  extend 
upstream  of  the  leading  edge.  The  leading-edge  separation  bubble  trapped  and 
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accumulated  seed  particles,  due  to  injection  from  the  reversed  flow.  This  seed  within  the 
bubble  was  able  to  sustain  itself  even  after  the  fog  was  turned  off.  The  bubble  also 
appeared  to  pulsate  in  size  in  accordance  with  the  intermittent  reversed-flow.  There  was 
some  three-dimensionality  to  the  reversed  flow  evidenced  by  small  spanwise  velocity 
traces.  Since  the  laser  sheet  was  only  projected  between  two  blades,  it  was  not  possible  to 
determine  whether  the  cascade  was  experiencing  intermittent  stall  from  blade  passage  to 
blade  passage,  or  if  the  entire  cascade  was  stalling  simultaneously. 


37 


F.  COMPUTATIONAL  RESULTS 


1.  Comparison  of  Blade-Surface  Pressure  Distributions 

The  predicted  and  experimental  blade-surface  pressure  distributions  are  shown  for 
40,  43,  46,  and  48  degrees  in  Figures  24-27  respectively.  Experimental  data  for  the  40 
and  43  degree  case  were  from  Reference  14.  The  46  degree  case  used  Reference  15 
experimental  data  and  the  48  degree  case  used  Reference  3  data.  The  AVDR  parameter 
was  held  constant  at  6%  for  all  solutions,  which  showed  good  agreement  for  the  first  three 
cases.  The  48  degree  prediction  over-estimated  the  experimental  values  aft  of  60%  chord. 
This  was  most  probably  due  to  under-estimating  the  stream-tube  thickness  variation  at  this 
inlet  flow  angle.  Increasing  the  AVDR,  as  inlet  flow  angle  increased,  improved  the 
agreement  of  the  computational  results;  however,  this  was  not  the  case  experimentally,  as 
documented  in  References  14,  16,  and  17. 

The  computed  static-pressure  recoveiy  far  downstream  of  the  trailing  edge  and 
that  measured  (from  Reference  18)  showed  good  agreement  for  the  40  and  43  degree 
cases.  The  computed  value  for  40  degrees  was  0.313  as  compared  with  a  measured  range 
of  0.32  -  0.34.  The  computed  value  for  43  degrees  was  0.357  as  compared  with  a 
measured  range  of  0.37  -  0.39.  The  computed  static-pressure  rise  far  downstream  is 
relevant  to  accuracy  in  predicting  the  losses  using  viscous  codes,  and  it  was  shown  here  to 
agree  well  with  experiments  near  design  incidence. 

For  Pi  =  48  degrees,  the  predicted  flow  field  pressure  contours  between  the  blade 
rows  are  shown  in  Figure  28.  The  pressure  field  was  continuous  at  the  grid  boundaries 
showing  good  periodicity;  however,  there  was  some  indication  of  grid  contamination  of 
the  pressure  field  at  the  trailing  edge  cut.  The  suction  peak  and  stagnation  point  were 
both  clearly  evident.  The  static-pressure  rise,  or  adverse-pressure  gradient,  towards  the 
trailing  edge  on  the  suction  side  was  also  apparent.  Investigation  of  the  velocity  field  in 
this  region  showed  that  there  was  flow  separation  aft  of  mid-chord,  followed  by  re¬ 
attachment  prior  to  the  blunt  trailing  edge.  Similar  behavior  occurred  for  the  46  degree 
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solution.  Experimentally  this  was  not  the  case,  as  the  flow  remained  attached  on  the  aft 
section  of  the  suction  surface.  This  indicated  that  the  code  may  be  predicting  incipient 
stall  at  48  degrees  (or  at  a  lower  inlet-flow  angle  of  46  degrees);  whereas,  experimentally 
this  was  only  found  at  50  degrees.  It  appeared  that  an  unsteady  solution  may  have  been 
computed  as  evidenced  by  the  climbing  density  residuals,  which  are  shown  in  Appendix  G, 
and  the  coalescence  of  the  pressure  contour  lines  throughout  the  flow  field  (Figure  28). 
The  40  and  43  degree  solutions  did  not  show  any  suction-side  separation  on  the  aft 
portion  of  the  blade,  in  agreement  with  experiments. 
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Figure  24.  40  Degree  Pressure  Distribution 
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Figure  26.  46  Degree  Pressure  Distribution 


T 


T 


^  -p  ^  ^ 


Cray  48.29  degrees  — 
Exp:  48  Press  surf  o 
Exp:  48  Suction  surf  + 


1 


I 


1 


0.2  0.4  0.6  0.8 

X/C 


Figure  27.  48  Degree  Pressure  Distribution 
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2.  Leading-Edge  Separation  Bubble 

A  leading-edge  separation  bubble  was  predicted  for  all  computational  solutions. 
The  location  of  the  computed  separation  point  was  insensitive  to  varying  inlet  flow  angle. 
The  re-attachment  point  for  each  case  was  determined  graphically.  Figure  29  shows  the 
bubble  region  for  Pi  =48  degrees  in  which  the  recirculating  flow  was  evident.  A 
comparison  of  computed  and  experimental  re-attachment  points  is  shown  in  Figure  30. 
The  experimental  data  was  from  Reference  19,  obtained  using  the  china-clay  technique. 
The  plot  showed  an  under-prediction  of  greater  than  50%  for  the  re-attachment  locations, 
which  is  typical  of  bubble  computations.  The  Baldwin-Lomax  turbulence  model  does  not 
account  for  freestream  turbulence,  which,  as  pointed  out  in  Reference  3,  is  necessary  in 
order  to  properly  characterize  the  development  of  the  blade  surface  boundary  layer,  or  the 
occurrence  of  separation.  The  two  most  likely  reasons  for  the  shortcoming  of  the 
computations  are  firstly,  the  limitations  of  the  Baldwin-Lomax  turbylence  model  and 
secondly,  the  lack  of  grid  refinement  in  the  streamwise  direction  aft  of  the  leading  edge. 
The  trend  of  increasing  downstream  distance  for  re-attachment,  as  inlet  flow  angle 
increased,  appeared  to  be  properly  computed. 
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Figure  29.  Leading-Edge  Velocity  Plot 
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3. 


Losses 


The  final  assessment  of  RVCQ3D  was  its  loss-prediction  capability.  Figure  3 1 
plots  measured  losses  as  reported  in  References  14,  16  and  17  with  those  computed  in  the 
present  study.  The  first  two  computed  cases  of  40  and  43  degrees  were  within  5%  of  the 
measured  results.  The  46  degree  case  showed  a  sharp  departure  from  the  measured 
results,  falling  short  by  45%.  An  attempt  to  increase  the  AVDR  parameter  for  this  case 
led  to  a  small  increase  in  the  loss  prediction  but  with  a  marked  degradation  of  the  blade 
pressure  distribution. 

The  output  file  for  each  solution  contained  a  listing  of  the  Baldwin-Lomax 
turbulence-model  behavior  by  printing  values  of  turbulent  eddy  viscosity  for  each  grid 
point.  Upon  interrogation  of  this  output  file,  the  behavior  of  the  Baldwin-Lomax 
turbulence  model  was  seen  to  exhibit  certain  characteristics.  Firstly,  for  all  cases  the 
turbulence  model  turned  on  inside  the  separation  bubble  region.  Secondly,  for  the  40 
degree  case,  the  model  did  not  turn  on  until  two  grid  points  inside  the  bubble  region, 
indicating  a  laminar  region  in  the  forward  section  of  the  bubble.  The  model  then  turned 
off  again  at  the  re-attachment  point  and  turned  back  on  further  downstream.  Thirdly,  for 
48  degrees  the  model  turned  on  at  the  separation  point  and  remained  on  past  the  re¬ 
attachment  point  to  the  trailing  edge.  The  transition  criteria  used  for  all  the  computations 
was  that  once  the  turbulence  model  computed  a  turbulent  eddy  viscosity  (|j,t)  which  was 
14  times  larger  than  the  laminar  viscosity  (hl),  this  value  of  pr  was  used  to  compute  an 
effective  viscosity  (peff  =  Pl  +  tix  )•  Otherwise  the  effective  viscosity  was  simply  the 
laminar  viscosity. 
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Figure  31.  Loss  Comparison 


47 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Experimental  flow  field  measurements  near  stall  were  successfully  repeated  for  a 
controlled-difiusion  compressor  cascade  using  a  two-dimensional  laser-Doppler 
velocimeter  system.  The  nominal  inlet  flow  angle  achieved  was  50.05  degrees  at  a 
constant  Mach  number  of  0.25  and  Reynolds  number  of  700,000.  The  measurements 
were  useful  in  that  they  showed  repeatability  with  previous  work  by  Ganaim,  and 
validated  a  test  case  for  two-dimensional  unsteady  separated  flows. 

Analysis  of  the  data  was  presented  through  a  unique  editing  process  which 
captured  the  character  of  the  intermittent  reversed  flow.  This  process  more  properly 
represented  the  stalling  region  as  a  combination  of  primarily  positive  and  negative  flow, 
since  the  time-mean  velocity  rarely  occurred  in  the  intermittent  region. 

Power-spectrum  and  auto-correlation  calculations  were  completed  on  the  data  for 
eleven  locations.  Interpretation  of  these  results  was  limited  to  their  general  behavior  due 
to  averaging  of  the  data  by  the  software  package.  Three  locations  were  chosen  for 
analysis  and  the  results  were  found  to  depend  greatly  on  the  data  acquisition  rate.  When 
expecting  to  perform  power-spectrum  analysis  on  experimental  data,  it  is  important  to 
ensure  that  the  data  rate  is  high  by  repeated  cleaning  and  proper  re-fitting  of  the  optical- 
access  window. 

Numerical  computations  showed  good  predictability  of  the  blade  pressure 
distribution  up  to  46  degrees  (approximately  6  degrees  above  design  incidence).  The  48- 
degree  case  began  to  deviate  by  up  to  50%  fi’om  measured  values  aft  of  two-thirds  chord. 
A  case-by-case  basis  for  adjusting  AVDR  may  be  necessary  to  bring  the  computed  results 
into  better  agreement  with  the  measurements  at  higher  inlet  flow  angles. 

A  leading-edge  separation  bubble  was  formed  in  all  computations  in  accordance 
with  experiments.  The  extent  of  the  computed  bubble  fell  short  of  that  seen  in 
measurements.  Computed  losses  agreed  well  at  low  incidence  angles.  As  the  inlet  flow 
angle  increased  above  46  degrees,  predicted  losses  deviated  by  45%  from  the 
measurements.  The  use  of  a  constant  AVDR  for  numerical  solutions  greatly  simplified  the 
process.  This  was  found  to  be  acceptable  for  analyzing  near  on-design  inlet  flow  angles. 
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Future  computations  of  the  flow  through  the  cascade  with  a  quasi-3d  code  should 
focus  on  achieving  the  following; 

1.  Continue  numerical  solutions  to  determine  the  optimum  AVDR  as  inlet-flow 
angle  varies,  and  compare  with  those  measured  in  the  cascade. 

2.  Quantify  what  effects  adjusting  AVDR  to  match  pressure  distributions  has  on 
the  computed  losses. 

3.  Increase  the  grid  size  in  both  the  I  and  J  directions  for  better  characterization  of 
the  flow  near  the  blade  surface  and  around  the  leading  edge  separation  bubble. 

4.  Replace  the  Baldwin-Lomax  turbulence  model  with  a  two-equation  turbulence 
model  to  account  for  fireestream  turbulence  intensity. 

5.  Perform  time-accurate  solutions  at  higher  incidence  angles  to  see  if 
unsteadiness  is  being  properly  modeled. 
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APPENDIX  A;  EXPERIMENTAL-MEASUREMENTS  LOG 


A  detailed  log  of  the  atmospheric  conditions,  tunnel  plenum  chamber 
conditions  and  LDV  traverse-mechanism  geometry  for  each  experimental  run. 


APPENDIX  B:  CALCULATION  OF  INLET-FLOW  REFERENCE  VELOCITY 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


OUTPUT  FROM  PROGRAM  CALIBRATE 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

LEAST  SQUARES  STRAIGHT  LINE  CURVE  FIT  IS  USED 
TO  DETERMINE  TUNNEL  CHARACTERISTICS  AT  DIFFERENT  SPEEDS 

NEWTON  S  METHOD  IS  USED  TO  DETERMINE  THE  REFERENCE  VELOCITY 
FROM  THE  RECORDED  AMBIENT  PRESSURE  AND  TUNNEL  PLENUM 
PRESSURE  AND  TEMPERATURE 

BEGIN  DETERMINING  TUNNEL  CHARACTERISTICS 
FROM  THE  FOLLOWING  MEASURED  VALUES 

AXIAL  VEL  TANGENTIAL  VEL  AMBIENT  PRESS.  PLENUM  PRESS.  PLENUM  TEMP. 

M  PER  SEC.  M  PER  SEC.  INCHES  MERCURY  INCHES  WATER  DEG.  C. 


18.3100 

22.1350 

29.9702 

1.7000 

13.8888 

32.4070 

38.5280 

29.9702 

4.7000 

13.8888 

42.5440 

50.1820 

29.9702 

8.1000 

15.5555 

47.2920 

56.3530 

29.9702 

10.0000 

16.6666 

51.7870 

61.9370 

29.9702 

12.1000 

17.2222 

55.3660 

66.7010 

29.9702 

14.0000 

18.3333 

CALCULATED  VALUES  FOR  THE  TUNNEL  CONFIGURATION 


TOTAL  VELOCITY  MACH  NUMBER 


MACH  NUMBER  FUNCT.  PRESSURE  RATIO 


0.2872654391 OE+02 
0.5034501 3984E+02 
0.65789247298E+02 
0.73567614295E+02 
0.80734660079E+02 
0.866857-39064E+02 


0.3781 8777200E-01 
0.66279705382E-01 
0.86361830486E-01 
0.963872351 72E-01 
0.1 056761 5343E+00 
0.11 3249251 18E+00 


0.49880294974E-02 
0.15207192279E-01 
0.25620261319E-01 
0.317667591 19E-01 
0.38003966212E-01 
0.43463396462E-01 


0.43443528269E-02 

0.11919477085E-01 

0.20366465634E-01 

0.25024235973E-01 

0.30121036781E-01 

0.34686726827E-01 


CALLING  LEAST  SQUARES  SUBROUTINE 
TO  DETERMINE  THE  PRESSURE  RATIO  AS  A  FUNCTION  OF  MACH  NO.  PARAM 


PRESSURE  RATIO  =  A1  *  ANUX  +  AO 


MATRIX  EQUATION 

0.60E+01  0.16E+OOAO  0.13E+00 
0.16E+00  0.53E-02A1  0.42E-02 

A1  =  0.78899597626E+00  A0  =  0.16213280706E-03 
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REFERENCE  CONDITIONS  FOR  EACH  RUN 


AMB  PRESS 
INCHES  Hg 

PLENUM  PRESS 
INCHES  H20 

PLENUM  TEMP 
DEGREES  C 

RUN  NAME 

29.9295 

11.9000 

18.3333 

1114s1  A.PRN 

30.1331 

12.0000 

14.4444 

ST1  3K.PRN 

30.1331 

12.0500 

14.4444 

ST1  30K.PRN 

30.1331 

12.0000 

15.5555 

ST1  50K.PRN 

30.1127 

12.0500 

16.6666 

1121S  1B.PRN 

30.1127 

12.1000 

18.3333 

1121S  1D.PRN 

30.2349 

11.9500 

16.9444 

1122S  2A.PRN 

30.1331 

12.0000 

17.7777 

ST2A  3K.PRN 

29.9295 

12.1500 

14.4444 

ST2A  30K.PRN 

30.1942 

11.8500 

17.2222 

1123S  3.PRN 

30.1331 

12.0000 

16.6666 

ST37  3K.PRN 

29.9295 

12.0000 

16.6666 

ST37  30K.PRN 

30.2145 

12.0000 

20.0000 

1123S  7.PRN 

30.1331 

12.1000 

15.5555 

ST73  3K.PRN 

29.9295 

12.0000 

17.2222 

ST73  30K.PRN 

30.1331 

12.1000 

14.4444 

ST7  3K.PRN 

30.0313 

11.9000 

14.4444 

ST7  30K.PRN 

30.0517 

12.1000 

17.2222 

ST10  EDG.PRN 

30.0313 

12.0500 

16.6666 

ST10  OUT.PRN 

30.1738 

12.1000 

15.6566 

1127S  14.PRN 

30.2145 

12.0500 

16.1111 

1128S  15.PRN 

30.2145 

12.1000 

17.2222 

1128S  17.PRN 

30.0313 

11.9000 

16.1111 

ST17  IN.PRN 

30.0313 

11.9000 

16.6666 

ST17  OUT.PRN 

30.2349 

11.9000 

18.3333 

1128S  18.PRN 

30.2145 

12.0000 

19.4444 

1128S  19.PRN 

30.0517 

12.2500 

16.6666 

ST19  IN.PRN 

30.0517 

1=  1 

12.1000 

16.6666 

ST19  OUT.PRN 

PRESSURE  RATIO  =  0.02968  MACH  NUMBER  PARAMETER  =  0.3741E-01 
RUN  NAME  =  1114s1  A.PRN 


BEGIN  NEWTON  ITERATION 


ITERATION  NUMBER 
ITERATION  NUMBER 
ITERATION  NUMBER 
ITERATION  NUMBER 


1  MACH  NO.  PARAM. 

2  MACH  NO.  PARAM. 

3  MACH  NO.  PARAM. 

4  MACH  NO.  PARAM. 


0.105475  ERROR  TERM 
0.104822  ERROR  TERM 
0.104821  ERROR  TERM 
0.104821  ERROR  TERM 


0.6522E-03 

0.1641E-05 

4).3698E-09 

0.8591E-13 


VREF=  80.23415022742 


1=  28 

PRESSURE  RATIO  =  0.03004  MACH  NUMBER  PARAMETER  =  0.3787E-01 
RUN  NAME  =  ST19  OUT.PRN 


BEGIN  NEWTON  ITERATION 


ITERATION  NUMBER  1  MACH  NO.  PARAM. 
ITERATION  NUMBER  2  MACH  NO.  PARAM. 
ITERATION  NUMBER  3  MACH  NO.  PARAM. 


0.105777  ERROR  TERM 
0.105485  ERROR  TERM 
0.105485  ERROR  TERM 


0.2921  E-03 
0.2872E-06 
-0.681  OE-10 


VREF=  80.51152210108 


54 


FILE  NAME 


EXPERIMENT 

NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 


REFERENCE 

VELOCITY 

80.2342 

79.7607 

79.9273 

79.9146 

80.2627 

80.6606 

79.8038 

80.2216 

80.5328 

79.5598 

80.0682 

80.3411 

80.4184 

80.2481 

80.4180 

80.0935 

79.5615 

80.5887 

80.3718 

80.1937 

80.0501 

80.3705 

79.7917 

79.8682 

79.8265 

80.3422 

81.0107 

80.5115 


1114s1_A.PRN 
ST1_3K.PRN 
ST1_30K.PRN 
ST1_50K.PRN 
1121S  1B.PRN 
1121S_1D.PRN 
1122S_2A.PRN 
ST2a3k.PRN 
ST2A  30K.PRN 
1123S_3.PRN 
ST37  3K.PRN 
ST37_30K.PRN 
1123S_7.PRN 
ST73_3K.PRN 
ST73_30K.PRN 
ST7_3K.PRN 
ST7_30K.PRN 
ST10_EDG.PRN 
ST10_OUT.PRN 
1127S_14.PRN 
1128S_15.PRN 
1128S_17.PRN 
ST17JN.PRN 
ST17_OUT.PRN 
1128S_18.PRN 
1128S_19.PRN 
ST19_IN.PRN 
ST19  OUT.PRN 
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APPENDIX  C:  VELOCITY  HISTOGRAMS,  TIME-HISTORY  PLOTS,  POWER- 
SPECTRUM  PLOTS,  AND  AUTO-CORRELATION  PLOTS 

1.  Station  3,  Point  7 

a.  3K  Data  Points 
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Uelocity  (w/s)  J  Uelocity  (n/s) 


Pos.  :  (-0.8394,  0.0001,  -4.5428)  st37_3k.i 

V  Component 

112.972  -1 - - - — I - 


-80.912  -H - 

O.OOOE«00 


89.355 


-62.149  H - 

O.OOQE^OO 


Component 


s.  :  (-0.8394,  0.0001,  -4.5428)  st37_3k. 

V  Component 

112.972  -1 - — - 1 - 


-88.912  H - 

O.OOOE^OO 


89.355 


-62.149  H - 

O.OOOE+00 


U  Component 


POWER  SPECTRUM  POWER  SPECTRUM 


CORRELATION  CORRELATION 


b.  30K  Data  Points 


Filename:  c:Sldw50a jwSst3Sst37_30k.s01  #  Processors:  2 

10.000  A - * - ' - 1 - ' - 


7.500 


Mode:  Coincidence 


-242.003 


-laoioiO  -35;S33  67.552  170-737 

Uelocity  (m/s)  for  Component  M 


Position  (rec/in)  ("”0-0393#  0.0000#  4-5420) 

Velocity  mean  =  15.972 


Velocity  at  cursor 
Percent  at  cursor  = 


273. y 


0.2 

2.64 


Filename:  c:Mdw50ajwSst3Sst37_30k.s01  »  Processors:  2 

10.000  A - ' - ' - 1 - - 


Mode •  Co i nc idence 


7.500 


5.000 


2.500 


0-000  -4 - 

-310-658 


-178.881  -47.185  84.672  216.449 

Velocity  (m/s)  for  Component  V 


Position  (rec/in)  (—0-8393#  0.0000#  4.5420) 

Velocity  mean  =  18.801 


Velocity  at  cursor 
Percent  at  cursor  = 


348.2 


0.1 

1.80 


Pos .  ' 
7  117, 


(-0.8393, 


0-0000, 

V 


-4.5420)  st37. 
Component 


301c.  vOl 


628 


-88.052 


O.OOOE^OO 


ms 

M  Compofient 


1.939E+05 


93.371 


^  -61 


u 

o 


o 


434 
O.OOOE+00 


Edit  Zoom  Hel 


Total  Tine  (ms)  =  60012 
Data  Rate  (KHz)  =  0.045 
Hean  Uelocity  =  18.78 
Standard  Deu.  =  32.94 
Uel.  at  Cursor  =  25.38 
Delta  Time  (ms)  =  700.8 
Tine  at  Cursor  =  10150 


Data  Rate  (KHz)  =  8.045 
Hean  Uelocity  =  15.96 
Standard  Dev.  =  25.88 
Uel,  at  Cursor  =  19.03 


1.939E+85 


Pos.  '  (-0.8393, 


0.0000, 


7  117.628 


-88.052 


-4.5428)  st37. 
Component 


OOk.vOl 


-TTt  Ifl/IJ 

Twj  Ml  uLlln/- 

O.OQOE+OO 


93.371 


ms 

M  Component 


1 .385E+83 


$  -61.434 


u 

o 


— tt 

n  mnfin  7 

f  1  (  [  LI  U  LAU 

I  Edit  Zoom  Hel 

Total  Time  (ns)  =  60012 
Data  Rate  (KHz)  =  0.845 
Hean  Uelocity  =  18.78 
Standard  Dev.  =  32.94 
Uel.  at  Cursor  =  73.08 
Delta  Tine  (ms)  =  5.008 
Time  at  Cursor  =  725.8 


Data  Rate  (KHz)  =  0.845 
Hean  Uelocity  =  15.96 
Standard  Dev.  =25.88 
Uel.  at  Cursor  =  55.25 


-O.OOOE+OO 


ms 


1 .385E+03 
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POWER  SPECTRUM  POWER  SPECTRUM 


File:c:Siau5OajwVst3Sst37_3OR.M01  TitlerPower  Spectrum 


10 

FREQUEHCV  (Hz) 


Graph  Component •  U 

y.  of  Bins  Containing  Data  =  100 

Scaling  Factors  (XjV)  =  1.0  ,  1.0 


File:c:Siaw50ajwSst3Sst37_30h.U01  TitlerPower  Spectrum 


FREQUENCy  (Hz) 

Graph  Component:  V 

y.  of  Bins  Containing  Data  =  100 

Scaling  Factors  (X^V)  =  1.0  ,  1.0 


CORRELATION  CORRELATION 


Station  7,  Point  3 


a.  3K  Data  Points 


Filename:  c :Sldu50a juSst7Sst73_3k .sOl  #  Processors:  2 
10.000  -I - - - > - 1 - ^ - * - 


7.500 


5.000 


2.500 


Node :  Co i nc idence 


-97.131  -52.396  -7.662  37.072  81.806 

Velocity  (m/s)  for  Component  14 


Position  (rec/in)  (-0.1225#  0.0000# 
Velocity  mean  =  14.705 


-3.2921)  Velocity  at  cursor 
Percent  at  cursor 


F  i  lename :  c :  S  ldu50a  jwSst7Sst73_3k .  sOl  It  Processors : 
10.000  - 1 - * - 1 - *- 


7.500 


Node:  Coincidence 


5.000 


2.500 


-190.119  -98.121  -6.124  85.874  177.872 

Velocity  (m/s)  for  Component  V 


Position  (rec/in)  (-0.1225#  0.0000#  -3.2921) 
Velocity  mean  =  39.875 


Velocity  at  cursor 
Percent  at  cursor  = 


Uelocity  («/s)  J  Velocity  (m/s) 


U  ComponcTit 


Total  Tiwe  (ms)  =  37455 
Data  Rate  (KHz)  =  0.082 
Mean  Velocity  =  39.87 
Standard  Deu.  =  22.99 
Mel.  at  Cursor  =  42.06 
Delta  Time  (ms)  =  55. OG 
Time  at  Cursor  =  7975. 


Data  Rate  (KHz)  =  0.002 
Hean  Velocity  =  14. 7G 
Standard  Dev.  =  11.18 
Vel.  at  Cursor  =  16.44 


L 


j.  :  (-0.1225,  0.0000,  -3.2921)  st73_3k.M01 

V  CompoTient 

100.873 


-29.123 

Q.OQQE^^OO 


laanmiiRai 


Total  Time  (ms)  =  37455 


Data  Rate  (KHz)  =  0.082 
Hean  Velocity  =  39.87 

Standard  Dev.  =  22.99 

Vel.  at  Cursor  =  37.21 
Delta  Time  (ms)  =  5.00Q 
1  3fl5E+03  Time  at  Cursor  =  725.0 


Data  Rate  (KHz)  =  0.082 
Hean  Velocity  =  14. 7G 
Standard  Dev.  =  11.18 
Vel.  at  Cursor  =  25.64 
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FOUER  SPECTRUM  POUER  SPECTRUM 


Fi le :c :\ldv50a ju\st7Sst73__3k .M01  Title :Power  Spectrum 


FREQUENCV  (Hz) 

Graph  Component :  U 
Z  of  Bins  Containing  Data  =  100 
Scaling  Factors  (X,Y)  =  1.0  ,  1.0 


^  File :c :Mdv50a jwSst7Sst73_3h . UOl  Title :Pouer  Spectrum 


Graph  Component:  V 

y,  of  Bins  Containing  Data  =  100 

Scaling  Factors  (XjY)  =  1.0  $  1.0 


FBEQUENCY  (Hz) 


CORRELATION  CORRELATION 


y,  of  Bins  Containing  Data  =  100 
Scaling  Factors  (X,Y)  =  1.0  f  1.0 


Scaling  Factors  (XjY)  =  1.0  f  1.0 
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b.  30K  Data  Points 


Filename:  c:Sldv50ajMSst7Sst73_301<.s01  #  Processors:  2 

10.000  -I - • - ' - 1 - ' - 


Hode2  Coincidence 


-^103.824 


-56,709  -9,595  37.519  84,633 

Uelocity  Cm/s)  for  Component  M 


Position  Crec/in)  C— 0.1224#  0.0001#  3.2920) 

Uelocity  mean  =  13,968 


Uelocity  at  cursor 
Percent  at  cursor  = 


131.7 


19.1 

1.45 


Filename:  c:Mdv5DajMSst7'sst73_301<.s01,#  Processors:  2  Mode:  Coincidence 

10.000  H - -• - ' - 1 -  I 


7.500 


5.000 


2.500 


0.000  - 

-200.597 


-110.138  -11.679  86.780  185.240 

Uelocity  (m/s)  for  Component  v 


Position  (rec/in)  (-0,1224#  0.0001#  3.2920) 

Uelocity  mean  =  37.561 


Uelocity  at  cursor 
Percent  at  cursor  = 


283.6 


48,3 

2.40 


Velocity  (m/s)  ^  J  Uelocity  (m/s) 


s.  :  (-0.1224,  0.0001,  -3.2920)  st73_30k.wOl 

V  Component 

111.394  -J - — 1 - 


-36.282  H - 

O.OOOE+00 


49.299 


-21.368  H - 

O.OOOE^OO 


U  Component 


-  Total  Time  (ms)  =  24086 

Data  Date  (KHz)  =  0.128 
Mean  Velocity  =  37.55 
Standard  Deu.  =  24.61 
Mel.  at  Cursor  =  35.04 

-  Delta  Time  (ms)  =  800.0 

2.216E+85  Time  at  Cursor  =  11600 


Data  Rate  (KHz)  =  0.128 
Mean  Velocity  =  13.96 
Standard  Dev.  =  11-77 
Vel,  at  Cursor  =  13.57 


2.216E+05 


5.  :  (-0.1224,  0.8001,  -3.2920)  st73_30k.v01 

V  Component 

111.394  -\ - —I - 


-36.282  H - 

O.O0OE+88 


49.299 


-21.368  H - 

O.OOOE+00 


M  Component 


-  Total  Time  (ms)  =  24086 

Data  Rate  (KHz)  =0.128 
I  Mean  Velocity  =  37.55 

M  ]  Standard  Dev.  =  24.61 

l|  ^  Vel.  at  Cursor  =  45.55 

-  Delta  Time  (ms)  =  5.000 

1.385E+03  Time  at  Cursor  =  725.0 


Data  Rate  (KHz)  =  0.128 
Mean  Velocity  =  13.96 
Standard  Dev.  =  11.77 
Vel,  at  Cursor  =  19.61 


1.385E+03 
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POWER  SFECTRUn  POWER  SPECTRUM 


Graph  Componept:  ^ 

X  of  Bins  Containing  Data  =  100 
Scaling  Factors  (X,V)  =  1.0  ,  1.0 


Graph  Component:  V 

Vi  of  Bins  Containing  Data  =  100 

Scaling  Factors  (X,V)  =  1.0  i  1-0 
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CORRELf^TION  CORRELATION 


1.000 


Fi le :c rSldvSOa jwSst7Sst73_30I< .UOl  Title rCorrelat ion  Graph 


yt  of  Bins  Containing  Data  =  100 
Scaling  Factors  (X,y)  =  1.0  i  1.0 


Scaling  Factors  (X,y)  =  1-0  ,  1.0 
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Ganaim  3K  Data  Points 


Filename:  c :Sldu__bert\st7S0507s71  .s03  #  Processors:  2 

10.000  - 1 - 1 - « - * - 


7.500 


Hode:  Coincidence 


5.000 


2.500 


-107.836  -58.575  -9.314  39.946  89.207 

Melocity  (m/s)  for  Component  U 


Position  (rec/in)  (-0.1224,  0.0000,  -3.2921) 
Melocitu  mean  =  15.316 


Uelocity  at  cursor 
Percent  at  cursor 


Filename:  c :Sldv_bertSst7S0507s71 .s03  It  Processors:  2 

10.000  H - * - * - 1 - * - 


7.500 


Mode:  Coincidence 


2  5.800 

£ 


2.500 


-222.621  -118.135  -13.648  98.839  195.325 

Velocity  (m/s)  for  Component  V 


Position  (rec/in)  (-0.1224,  0.0000,  -3.2921) 
Velocity  mean  =  38.608 


Velocity  at  cursor 
Percent  at  cursor 
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Pos.  :  (-0.1224,  0.0000,  -3.2921)  0507s71.v03 

V  Component 

116.950 


-39.742 


O.OOOE+00 


155E+03 


52.250 


W 

S  -21.633 


U 

o 

f-4 

V 


O.OOOE^OO 


4.155E+03 


I  Edit  Zoom  Hel 

Total  Time  (ns)  =  5064. 
Data  Rate  (KHz)  ==  0.606 
Hean  Uelocity  =  38.60 
Standard  Deu.  =  26.11 
Mel.  at  Cursor  =  0.113 
Delta  Tine  (ns)  =  15.00 
Tine  at  Cursor  =  2175. 


Data  Rate  (KHz)  =  0.607 
Mean  Uelocity  =  15.31 
Standard  Deu.  =  12.31 
Ue 1 .  at  Cursor  =  -2.67 


Pos.  :  (-0.1224,  0.0000,  -3.2921)  0507s71.v03 

V  Component 

116.958 


IfcdTir 


-39.742 


O.OOOE+00 
52.250 


ms 

U  (kimponent 


1.3052+03 


e 

E 


^  -21.633 


J 


J 


■wn 


O 

O 

0 

:> 


O.OOOE+00 


ms 


1.385E+03 


I  Edit  Zoom  Hel 

Total  Tine  (ns)  =  5064. 
Data  Rate  (KHz)  =  0.606 
Mean  Uelocity  =  38.60 
Standard  Dev.  =  26.11 
Uel.  at  Cursor  =  53.25 
Delta  Tine  (ns)  =  5.000 
Tine  at  Cursor  =  725.0 


Data  Rate  (KHz)  =  0.607 
Mean  Uelocity  =  15.31 
Standard  Dev.  =  12.31 
Uel.  at  Cursor  =  20.37 
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POWER  SFECTRUn  POWER  SPECTRUM 


CORRELATION  CORRELATION 


3, 


Station  7,  Point  20 


a.  3K  Data  Points 


Filename :  c :Mdw50a jMSst7\st7_3J< .sOl 

10.000  -J - 


Processors • 
J- 


Hode:  Coincidence 


c 

0) 

0 

u 

(S 


7.5G0  H 


5. GOO  H 


2. 500  H 


0.000 


-65.428 


-28!507  0.415  45.336  82.257 

Velocity  (m/s)  for  Component  U 


Position  (rec/in)  (0.3654,  0.0000,  3.2920) 

Velocity  mean  =  26.884 


Velocity  at  cursor 
Percent  at  cursor  = 


30.9 

3.22 


Filename:  c:Mdw50ajMSst7Sst7_3k.sOl  Processors!  2 

10.000  -I - • - — ^ - n 

7.500  - 

'P 

0) 

t  5.000  - 

e 

Oi 


Mode!  Coincidence 

- 1 - - 1- 


2.500  H 


0.000 


-165.506 


-75.660  14.i86  104.032  193.879 

Uelocitii  (m/s)  For  Component  V 


283.7 


Position  (rec/in)  (0.3654,  0.0000,  3.2920) 

Velocity  mean  =  59.130  _ 


Velocity  at  cursor  = 
Percent  at  cursor  = 


68.9 

1.73 
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Pos.  :  (0.3654,  O.OQOO,  -3.2920)  st7_3lc.u01 


O.OOOE+00 


54.541 


$  -O, 

0 

o 


773 
O.OOdE+00 


ms 

M  Component 


4.432E+04 


I  Edit  Zoom  Hel 

Total  Time  (ns)  =  44037 
Data  Rate  (KHz)  =  0.070 
Mean  Ueloclty  =  59.10 
Standard  Dev.  =  22.46 
Mel,  at  Cursor  =  68.00 
Delta  Time  (ms)  =  160.0 
Time  at  Cursor  23200 


Data  Rate  (KHz)  =  0.070 
Hean  Ueloclty  =  26.80 
Standard  Dev.  =  9.219 
Mel.  at  Cursor  =  29.90 


ms 


4.432E+04 


Pos.  :  (0.3654, 


126.494 


-8.275 


0.0000,  -3.2920)  st7_3k.v81 

M  Component 


O.OOOE+00 


U 


54.541 


ms 

Component 


1.385E^03 


(9 


•wi 

0 

o 

o 


-0.773 

O.OOOE+00 


71p3n; 

_ c-ws _ C 

-W- 

ri^ir 

1  "  ..y-Z] 

a 

Edit  Zoom  Hel 


Total  Time  (ns)  =  44037 
Data  Rate  (KHz)  =  0.070 
Hean  Melociiy  =  59.10 
Standard  Dev.  =  22.46 
Mel.  at  Cursor  =  74.07 
Delta  Time  (ms)  =  5.000 
Time  at  Cursor  =  725.0 


Data  Rate  (KHz)  =  0.070 
Hean  Melocity  =  26.80 
Standard  Dev,  =  9.219 
Mel.  at  Cursor  “  28.07 


ms 


1 .385E+83 


77 


FOUER  SPECTRUM  POWER  SPECTRUM 


F i I e : c : M  d  v5Ga  juSs t7Ss t7_3k • UDl  Title: Power  Spectrum 


FREQUENCV  (Hz) 

Graph  Component :  M 

V.  of  Bins  Containing  Data  =  100 

Scaling  Factors  (X,y)  =  1.0  j  1-0 


Fi le :c :Sldu50a jw\st7Sst7_3k -UOl  Title -Power  Spectrum 


FREQUENCy  (Hz) 

Graph  Component:  V 
X  of  Bins  Containing  Data  =  100 
Scaling  Factors  (X,y)  =  1.0  ,  1.0 


CORRELATION  CORRELATION 


1.000 


F  i  1  e :  c :  S 1  d  w50a  juSst7Sst7_31< .  UOl 


Title :Correlation  Graph 


Scaling  Factors  (X,y)  =  1.0  j  1.0 


Scaling  Factors  (X,Y)  =  1.0  ,  1.0 


79 


Filename;  c :Mdw50a jwSst7Nst7_30It .sOl  #  Processors:  2 

10.000  -I - ' - ' - 1 - ' - 


r.500 


Mode •  Co inc idence 


5.000 


2.500 


-201.488 


-99.572  2.345  104.261  206.178 

Uelocity  (m/s)  for  Component  V 


Position  (rec/in)  (0.3653j  — 0.8001 j  3.2928) 
Uelocity  mean  =  53.312 


Melocity  at  cursor 
Percent  at  cursor 


Uelociti^  (m/s)  ^  3*  Uelocity  (m/s) 


s.  :  (0.3653,  -0.0001,  -3.2920)  st7_30k.u01 

V  Component 

129.733  -A - - — I - 


-23.117  -I - 

D.000E*00 


-  Total  Time  (ms)  =  43243 

Data  Rate  (KHz)  =  0.071 

Mean  Velocity  =  53. 3G 

Z  Standard  Deu.  =  25.47 

Vel.  at  Cursor  =  60.62 

-  Delta  Time  (ms)  =  5.00G 

1.3S5E+03  Time  at  Cursor  =  725.0 


Data  Rate  (KHz)  =  0.071 
Mean  Uelocity  =  25.10 
Standard  Dev.  =  11.11 
Vel.  at  Cursor  =  31.40 
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POUER  SPECTRUn  POWER  SPECTRUM 


CORRELATION  CORRELATION 


Ganaim  3K  Data  Points 


Filename:  c :Sldv_bertSst7S0507s71 .s20  It  Processors:  2 

10.000  -\ - 1 - — - 1 - , - * - 


Mode:  Coincidence 


7.500 


5.000 


2.500 


0.000  -4 - 

~84 .334 


-39.979  4.376  48.731  93.086 

Ueloclty  (m/s)  for  Component  M 


Position  (rec/in)  (0.3654,  0.0000,  -3.2921) 
Uelocity  mean  =  26.554 


Me loci ty  at  cursor 
Percent  at  cursor  = 


137.4 


31.4 

4.82 


Filename:  c :Sldu_bertSst7\0507s71.s20  It  Processors:  2 

10.000  -A - ' - ' - 1 - • - 


7.500 


Mode:  Coincidence 


-208.803  -102.453  3.896  110,246  216.596 

Me  loci  ty  (m/s)  for  Compoinent  V 


Position  (rec/in)  (0.3654,  0.0000,  3.2921) 

Uelocity  mean  =  57,071 


Uelocity  at  cursor 
Percent  at  cursor  • 


ocity  (m/'s)  J  Uelocity  (m/s) 


5.  :  (0.3654,  0.0000,  -3.2921)  0507s71.u20 

V  Component 

136.033  H - - - — I - 


-22.691  H - 

O.OOOE+00 


- Total  Time  (ms)  =  5012. 

Data  Rate  (KHz)  =  0.613 

aj  Hean  Uelocity  =  57.07 

7B  standard  Dev.  =  26.58 

ft  Mel.  at  Cursor  =  39.96 

-  Delta  Time  (ms)  =  20.00 

5.540E+03  Time  at  Cursor  =  2900. 


Data  Rate  (KHz)  =  0.613 
Mean  Uelocity  =  26.55 
Standard  Dev.  =  11.08 
Mel.  at  Cursor  =  24.18 


(0.3654,  0.0000,  -3.2921)  0507s71.v20 


136 .033 


V  Component 


iiiiiiii«i!iiiiriiiiuru'iiiu«wvi[iiuii'~ . 


“22.691 

O.OOOE+00 


Edit  Zoom 


Total  Time  (ms)  =  5012. 
Data  Rate  (KHz)  =  0.613 
Mean  Uelocity  =  57.07 

Standard  Dev.  =  26.58 

Uel.  at  Cursor  =  77.31 
Delta  Time  (ms)  =  5.008 
1.385E+G3  Time  at  Cursor  =  725.8 


Data  Rate  (KHz)  =  0,613 
Mean  Uelocity  =  26.55, 
Standard  Dev,  =  11.08 
Uel.  at  Cursor  =  29.97 


POWER  SPECTRUM  POWER  SPECTRUM 


CORRELATION  CORRELATION 


Scaling  Factors  (X,Y)  =  1.0  #  1.0 


Scaling  Factors  (X#Y)  =  1.0  j  1.0 
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APPENDIX  D:  REDUCED  DATA  FROM  STATIONS  1,  3,  7, 19 


Pitchwise  Survey  at  Station  1 


X  (m)  I  y  m 


LOO 


-6.29 


U/Vr^ 


0.7560 


V/Vr^ 


0.6273 


UYMag/Vr^ 


0.9824 


U- Turbulence 


4.5167 


Yr  Turbulence 


4.2288 


l^teynpJd?^ 

Stress 


^.5645 


Cojrel^on 

CoeSici^ 


-0.0459 


0.75 


-6.29 


0.7527 


0.6348 


0.9847 


4.2846 


4.1545 


■0.1621 


■0.0141 


0.50 


-6.29 


0.7549 


0.6361 


0.9872 


4.1377 


4.1513 


■0.8599 


-0.0778 


0.25 


■6.29 


0.7532 


0.6393 


0.9879 


3.9633 


4.2387 


-0.0287 


-0.0027 


0.00 


-6.29 


0.7551 


0.6417 


0.9910 


4.2522 


4.0219 


-0.0638 


-0.0058 


-0.25 


-6.29 


0.7573 


0.6399 


0.9914 


4.1720 


4.0454 


-0.0756 


-0.0070 


-0.50 


-6.29 


0.7646 


0.6367 


0.9950 


4.3993 


3.9576 


0.4649 


0.0415 


-0.75 


-6.29 


0.7671 


0.6362 


0.9966 


4.3408 


3.9461 


0.1865 


0.0169 


-1.00 


-6.29 


0.7655 


0.6322 


0.9929 


4.5415 


3.8821 


0.4439 


0.0391 


-1.25 


-6.29 


0.7685 


0.6303 


0.9940 


4.0621 


3.6989 


0.3160 


0.0327 


-1.50 


-6.29 


0.7622 


0.6312 


0.9896 


4.9631 


3.7591 


0.0913 


0.0076 


■1.75 


-6.29 


0.7527 


0.6333 


0.9837 


5.6056 


3.5044 


0.2665 


0.0211 


-2.00 


-6.29 


0.7505 


0.6371 


0.9844 


4.0228 


3.4712 


-0.0511 


-0.0057 


-2.25 


-6.29 


0.7509 


0.6368 


0.9846 


3.6099 


3.5205 


-0.0113 


-0.0014 


-2.50 


-6.29 


0.7519 


0.6417 


0.9885 


3.5904 


3.3798 


0.1693 


0.0217 


-2.75 


-6.29 


0.7512 


0.6462 


0.9909 


3.5229 


3.4416 


-0.1665 


■0.0213 


■3.00 


-6.29 


0.7553 


0.6495 


0.9961 


3.5721 


3.2735 


0.1392 


0.0185 


■3.25 


-6.29 


0.7576 


0.6500 


0.9982 


3.9669 


3.3235 


-0.2241 


-0.0264 


-3.50 


-6.29 


0.7767 


0.6438 


1.0088 


3.5708 


3.4123 


-0.3948 


-0.0503 


-3.75 


■6.29 


0.7787 


0.6344 


1.0044 


3.5851 


3.5118 


-0.3790 


-0.0468 


-4.00 


-6.29 


0.7753 


0.6322 


1.0004 


3.4988 


3.3147 


-0.1981 


-0.0265 


-4.25 


-6.29 


0.7714 


0.6303 


0.9962 


3.4435 


3.3032 


-0.1314 


-0.0179 


50 

75 

00 

25 

50 

75 

00 

25 

50 

75 

00 

25 

50 

75 


-6.29 


0.7634 


0.6295 


0.9894 


3.7142 


3.1308 


0.0422 


0.0056 


-6.29 


0.7596 


0.6341 


0.9895 


3.4552 


2.9855 


0.2140 


0.0322 


-6.29 


0.7592 


0.6383 


0.9918 


3.3208 


2.9673 


-0.2577 


-0.0406 


-6.29 


0.7596 


0.6459 


0.9971 


2.9610 


2.8494 


-0.2546 


-0.0469 


-6.29 


0.7652 


0.6516 


1.0050 


3.2893 


2.7745 


0.0387 


0.0066 


-6.29 


0.7714 


0.6537 


1.0111 


2.9942 


3.0307 


-0.3735 


-0.0639 


-6.29 


0.7748 


0.6563 


1.0154 


3.2675 


2.9769 


-0.1074 


-0.0171 


-6.29 


0.7782 


0.6528 


1.0157 


2.9052 


2.6739 


-0.1011 


-0.0202 


-6.29 


0.7787 


0.6490 


1.0137 


3.4300 


2.7946 


-0.0564 


-0.0091 


-6.29 


0.7773 


0.6426 


1.0085 


2.7609 


2.5602 


-0.0539 


-0.0119 


-6.29 


0.7767 


0.6388 


1.0057 


3.0063 


2.5740 


0.4140 


0.0831 


-6.29 


0.7741 


0.6370 


1.0024 


3.1116 


2.4751 


0.2009 


0.0405 


-6.29 


0.7727 


0.6406 


1.0038 


2.7515 


2.3997 


0.2478 


0.0583 


-6.29 


0.7686 


0.6405 


1.0005 


2.5751 


2.2243 


0.1246 


0.0338 


-8.00 


■6.29 


0.7665 


0.6444 


1.0014 


2.8413 


2.2209 


0.1098 


0.0270 


-8.25 


-6.29 


0.7686 


0.6469 


1.0046 


3.1292 


2.2872 


0.0908 


0.0197 


-8.50 


-6.29 


0.7694 


0.6525 


1.0089 


3.0650 


2.0272 


0.0634 


0.0158 


-8.75 


-6.29 


0.7739 


0.6534 


1.0129 


2.9483 


2.1966 


0.0291 


0.0070 


-9.00 


-6.29 


0.7783 


0.6542 


1.0167 


2.4833 


2.1338 


0.0240 


0.0070 
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:  Total  Flow  Average 

Fo^ive  Average  | 

x{i») 

U/Vtef 

V/Vref 

Iftol/Vief 

U-Turt> 

V-Ttirb 

VtVnt 

V/Vief 

U/vtef 

V/Vref 

^.9160 

-4.5421 

-0.0508 

4).0820 

0.0964 

-6.8660 

-7.5121 

-0.0745 

-0.1005 

0.0976 

0.1242 

-0.9061 

-4.5420 

-0.0566 

-0.0813 

0.0991 

-7.0072 

-9.0429 

-0.0780 

-0.1061 

0.0952 

0.1348 

-0.8951 

-4.5420 

-0.0566 

-0.0807 

0.0986 

-7.3448 

-10.6218 

-0.0814 

-0.1091 

0.1169 

0.2125 

-0.8831 

-4.5420 

-0.0524 

-0.0713 

0.0885 

-10.4817 

-13.7893 

-0.0855 

-0.1150 

0.1897 

0.2506 

-0.8699 

-4.5420 

-0.0105 

-0.0184 

0.0212 

-18.3077 

-24.1420 

-0.0848 

-0.1154 

0.3524 

0.4773 

-0.8553 

-4.5420 

0.0692 

0.0767 

0.1033 

26.3741 

33.9905 

-0.0825 

-0.1146 

0.4479 

0.5715 

-0.8393 

-4.5420 

0.2345 

0.2896 

0.3726 

33.1398 

43.3500 

-0.0812 

-0.1178 

0.5142 

0.6610 

-0.8217 

-4.5420 

0.3414 

0.4153 

0.5377 

33.9309 

43.6112 

-0.0799 

-0.1182 

0.5525 

0.6905 

-0.8023 

-4.5420 

0.4373 

0.5191 

0.6787 

33.0556 

41.3025 

-0.0764 

-0.1118 

0.5858 

0.7099 

-0.7811 

-4.5420 

0.5376 

0.6040 

0.8086 

29.6847 

36.5104 

-0.0753 

-0.1088 

0.6291 

0.7162 

-0.7575 

-4.5420 

0.5532 

0.5886 

0.8078 

28.0982 

34.1066 

-0.0704 

-0.1176 

0.6245 

0.6757 

-0.7317 

-4.5420 

0.6236 

0.6658 

0.9122 

23.7893 

30.3195 

-0.0718 

-0.0978 

0.7168 

-0.7035 

-4.5420 

0.6598 

0.6660 

0.9375 

20.8838 

26.3703 

-0.0571 

-0.0992 

0.6841 

0.6953 

-0.6723 

-4.5420 

0.6963 

0.6885 

0.9792 

17.4786 

23.8027 

-0.0673 

-0.1251 

0.7080 

0.7041 

-0.6377 

-4.5420 

0.7192 

0.7057 

1.0076 

14.4658 

22.1652 

-0.0455 

-0.1565 

0.7255 

0.7169 

-0.6001 

-4.5420 

0.7387 

0.7207 

1.0320 

12.5289 

20.1208 

0.0000 

0.0000 

0.7422 

0.7275 

-0.5585 

-4.5420 

0.7553 

0.7316 

1.0515 

10.5465 

18.8519 

0.0000 

0.0000 

0.7565 

0.7352 

-0.5127 

-4.5420 

0.7781 

0.6664 

1.0245 

10.0197 

15.9950 

0.0000 

0.0000 

0.7794 

0.6687 

-0.4626 

-4.5420 

0.7642 

0.7640 

1.0806 

8.0755 

15.7110 

0.0000 

0.0000 

0.7642 

0.7643 

-0.4073 

-4.5420 

0.7645 

0.7642 

1.0809 

8.4760 

15,8142 

0.0000 

0.0000 

0.7645 

0.7647 

-0.3465 

-4.5420 

0.7524 

0.7933 

1.0933 

8.0195 

15.1845 

0.0000 

0.0000 

0.7524 

0.7933 

-0.2796 

-4.5420 

0.7442 

0.7777 

1.0764 

8.2561 

14.5299 

0.0000 

0.0000 

0.7777 

-0.2060 

-4.5420 

0.7460 

0.7736 

1.0746 

8.6669 

13.4328 

0.0000 

0.0000 

0.7460 

0.7736 

-0.1251 

-4,5420 

0.7330 

0.7705 

1.0635 

9.0455 

13.2214 

0.0000 

0.0000 

0.7330 

0.7705 

-0.0361 

-4,5420 

0.7322 

0.7560 

1.0525 

9.3525 

12.0596 

0.0000 

0.0000 

0.7322 

0.7560 

0.0619 

-4.5420 

0.7208 

0.7548 

1.0437 

9.5275 

11.5091 

0.0000 

0.0000 

0.7208 

0.7548 

0.1696 

-4.5420 

0.7186 

0.7403 

1.0317 

9.6073 

10.9092 

0.0000 

0.0000 

0.7186 

0.7406 

0.2882 

-4.5420 

0.6990 

0.7356 

1.0148 

9.5020 

10.1140 

0.0000 

0.0000 

0.6990 

i  0.7356 

0.4185 

-4.5420 

0.6903 

0.7245 

1.0007 

9.2803 

9.0998 

0.0000 

0.0000 

0.6903 

'  0.7245 

0,5620 

-4.5420 

0.6706 

0.7162 

0.9811 

9.0082 

8.4316 

0.0000 

0.0000 

0.6706 

0.7162 

0.7196 

-4.5420 

0.6449 

0.6946 

0.9479 

8.7474 

8.0618 

0.0000 

0.0000 

0.6449 

0.6946 

0.8932 

-4.5420 

0.6282 

1  0.6704 

0.9187 

7.8477 

7.6788 

0.0000 

0.0000 

0.6282 

0.6704 
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Total  Flow  Averaj^e 


Average  of  Negative  }  Average  of  Positive 


:  .  X  (my 

y/v^ 

y/vnef 

Utot/Vref 

U-Turb 

y-Tiwb 

U/Vrrf  .. 

v/v#ef 

,.v(yKt 

y/Vfef 

^.1435 

-3.2920 

0.1749 

0,4748 

0.5060 

12.2241 

26.7540 

4).0643 

-0.1126 

0.2064 

0.5479 

-0.1335 

-3.2920 

0.1809 

0.4820 

0.5148 

12.9712 

26.9304 

-0,0657 

-0.1081 

0.2146 

0.5524 

-0.1225 

-3.2920 

0.1780 

0.4843 

0.5159 

13.8963 

28.4924 

-0.0656 

-0.1100 

0.2211 

0.5748 

-0.1105 

-3.2920 

0.1861 

0.4919 

0,5259 

14.7300 

30.1213 

-0.0748 

-0.1079 

0.2333 

0.5898 

-0.0973 

-3.2920 

0.1889 

0.4992 

0.5338 

15.0800 

30.5573 

-0.0691 

-0.1077 

0.2385 

0.6005 

-0.0827 

-3.2920 

0.1857 

0.4874 

0.5216 

15.8231 

31.1735 

-0.0719 

-0.1014 

0.2387 

0.5935 

-0.0667 

-3.2920 

0.1938 

0.5026 

0.5387 

15.9031 

32.6549 

-0.0739 

-0.1024 

0.2487 

0.6161 

-0.0481 

-3.2920 

0.2026 

0.5130 

0.5515 

16.1680 

33.0578 

-0.0707 

-0.1090 

0.2588 

0.6283 

-0.0297 

-3.2920 

0.2116 

0.5290 

0.5698 

16.4695 

33.2063 

-0.0736 

-0.0988 

0.2664 

0.6396 

-0.0084 

-3.2920 

0.2166 

0.5227 

0.5658 

16.8288 

35.1634 

-0.0681 

-0.0991 

0,2776 

0.6503 

0.0151 

-3.2920 

0.2136 

0.5201 

0.5622 

16.7788 

36.0857 

-0.0675 

-0.0983 

0.2728 

0.6490 

0.0409 

-3.2920 

0.2162 

0.5202 

0.5634 

16.9932 

36.5286 

-0.0664 

-0.0981 

0.2778 

0.6529 

0.0692 

-3.2920 

0.2401 

0.5600 

0.6093 

17,0275 

37.3337 

-0.0654 

-0.0937 

0.3014 

0.6985 

0.1004 

-3.2920 

0.2453 

0.5712 

0.6216 

16.7892 

37.0429 

-0.0686 

-0.0993 

0.2974 

0.6905 

0.1348 

-3.2920 

0.2533 

0.5765 

0.6297 

16.9775 

38.0062 

-0.0604 

-0.0986 

0.3091 

0.7065 

0.1725 

-3.2920 

0.2580 

0.5909 

0.6448 

16.8549 

37.0968 

-0.0574 

-0.0974 

0.3106 

0.7097 

0.2141 

-3.2920 

0.2733 

0.6101 

0.6685 

16.1963 

36.6200 

-0.0642 

-0.0975 

0.3161 

0.7126 

0.2599 

-3.2920 

0.2849 

0.6250 

0.6869 

15.8988 

36.1124 

-0.0660 

-0.0930 

0.3271 

0.7290 

0.3100 

0.2988 

0.6631 

0.7273 

15.0613 

33.7251 

-0.0568 

-0.0940 

0.3324 

0.7409 

0.3654 

-3.2920 

0.3156 

0.7053 

0.7726 

14.0393 

31.8136 

-0.0652 

-0.0853 

0.3420 

0.7700 

0.4259 

-3,2920 

0.3282 

0.7261 

0.7969 

12.8077 

29.6575 

-0.0678 

-0.0934 

0.3481 

0.7733 

0.4930 

-3.2920 

0.3323 

0.7293 

0.8014 

12.6994 

27.7864 

-0.0715 

-0.1135 

0.3496 

0.7683 

0.5666 

-3.2920 

0.3469 

0.7369 

0.8145 

11.3452 

25.8805 

-0.0583 

-0.0838 

0.3550 

0.7616 

0.6474 

-3.2920 

0.3549 

0.7283 

0.8102 

11.2177 

24.6712 

-0.0594 

-0.0810 

0.3644 

0.7534 

0.7365 

-3.2920 

0.3667 

0.7457 

0.8310 

10.4527 

22.1508 

-0.0553 

-0.0619 

0.3716 

0.7615 

0.8345 

-3.2920 

0.3819 

0.7624 

0.8527 

9.8618 

18.9547 

0.0000 

0.0000 

0.3835 

0.7707 

0.9422 

-3.2920 

0.3868 

0.7551 

0.8484 

10.0925 

17.5833 

0.0000 

-0.0553 

0.3885 

0.7614 

1.0608 

-3.2920 

0.3912 

0.7689 

0.8627 

9.6454 

15.1138 

0.0000 

1  0.0000 

0.3914 

0.7730 

1.1911 

-3.2920 

0.3919 

0.7749 

0.8683 

9.1682 

12.5425 

0.0000 

;  0.0000 

0.3921 

0.7761 

1.3345 

-3,2920 

0.3849 

0.7814 

0.8710 

8.8570 

10.8529 

0.0000 

;  0.0000 

0.3849 

0.7814 

1.4922 

-3.2920 

0.3739 

!  0.7863 

0,8706 

7.9523 

8.9572 

0.0000 

0.0000 

0.3739 

0.7863 

1.6657 

-3.2920 

0.3500 

0.7726 

0.8482 

'  6.8941 

7.5859 

0.0000 

0.0000 

0.3500 

0.7729 
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X  (niches)  |Y  Cmc^es)  U/  Vref  V/Vief 


Average  of  Negative  |  Average  of  Positive 


U/Vtef  V/Vief  U/Vr^  V/Vref 
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APPENDIX  E:  CFD  GRID  GENERATION 


The  FORTRAN  grid  generation  program  employed  was  the  modified  version  of 
GRAPE  (GRids  About  Aerofoils  using  Poisson’s  Equation)  by  R.  V.  Chima,  [Ref  11]. 
Prior  to  this  study,  the  grid’s  leading  edge  was  blunt  or  flat-faced,  which  led  to 
inaccuracies  in  the  computed  solutions.  It  was  determined  that  GRAPE  properly 
smoothed  the  points  around  the  leading  edge,  and  that  it  required  definition  of  the  leading 
edge  with  more  user-supplied  coordinates.  Careful  attention  was  given  to  properly 
blending  the  top  and  bottom  grid  surfaces  into  the  0.045  inch  radius  circular  arc  that 
defines  the  leading  edge.  This  was  done  by  ensuring  that  first  derivative  continuity  was 
within  a  few  percent  difference.  The  final  coordinates  for  the  CD  blade  definition  are 
contained  below  in  namelist  &GRID3. 


&GRID1 

JMAX=240,KMAX-49,NTBTYP=3,NAIRF=5,JAIRF=163,NIBDST=?, 

DSI=0.0002,ITEBOT=60,JTETOP*181,XLE=OANOBSHP=7,XTE=1.0, 

ALAMF=O.G,ALAMR*0.0,XLEFT=-.3,XRIGHT=3., 

RCORN-  l0,NOUr=4,NORDA=5,3,>IAXITA»l00,100, 

&END 

&GR1D2 

AAAI=.3,BBBI=3,DSOBI=0.01,ROTANG=^14.4, 
XTFRAC=I.5,PITCH=.5988,DSRA=.4920,DSLE=000S,DSTE=.O0l, 
NLE=-2a4TO=i8,WAKEP=l.,OMEGR=1.0,OMBGS=l  .0,OMEGP=1 .0,OMEGQ=1 .0, 
&END 

&GRID3  AIRFX== 

0.9996,  0.9971,  0.9929,  0.9876,  0.9830,  0.9803,  0.9776,  0.9749, 

0.9601,  0.9453,  0.9305,  0.9157,  0.9010,  0.8862,  0.8714,  0.8567, 

0.8419,  0:8271,  0.8124,  0.7976,  0.7828,  0.7681,  0.7533,  0.7385, 

0.7238,  0.7090,  0.6942,  0.6795,  0.6647,  0.6499,  0.6352,  0.6204, 

0.6056,  0.5908,  0.5760,  0.5612,  0.5465,  0.5317,  0.5169,  0.5022, 

0.4874,  0.4726,  0.4579,  0.4431,  0.4283.  0.4136,  0.3988,  0.3840, 

0.3693,  0.3545,  0.3397,  0.3250,  0.3102,  0.2954,  0.2807,  0.2659, 

0.2511,  0.2364,  0.2216,  0.2068,  0.1920,  0.1772,  0.1624,  0.1477, 

0.1329,  0.1181,  0.1034,  0.0886,  0.0738,  0.0591,  0.0443,  0.0343, 

0,0244,  0.0114,  0,0097,  0.0080,  0.0063,  0.0047, 

0.0033,  0.0021,  0.0011,  0.0005,  0.0001,  O^OOOl,  0.0004, 

0.0010,  0.0019,  0.0030,  0.0044,  0.0045,  0.0177,  0.0310,  0.0443, 

0.0591,  0.0738,  0.0886,  0.1034,  0.1181,  0.1329,  0.1477,  0.1624, 

0.1772,  0.1920,  0.2068,  0.2216,  0.2364,  0.2511,  0.2659,  0.2807, 

0.2954,  0.3102,  0.3250,  0.3397,  0.3545,  0.3693,  0.3840,  0.3988, 

0.4136,  0.4283,  0.4431,  0.4579,  0.4726,  0.4874,  0.5022,  0.5169,  .  .. 

0.5317,  0.5465,  0.5612,  0.5760,  0.5908,  0.6056,  0.6204,  0.6352, 
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0.6499,  0.6647,  0.6795,  0.6942, 
0.7681,  0.7828,  0.7976,  0.8124, 
0.8862,  0.9010,  0.9157,  0.9305, 
0.9855,  0.9908,  0.9956,  0.9988, 
A1RFV= 

0.0093,  0.0044,  0,0012,  0.0000, 
0.0089,  0.0135,  0.0178,  0.0218, 
0.0381,  0.0407,  0.0431,  0.0453, 
0.0546,  0.0561,  0.0575,  0.0589, 
0.0643,  0.0652,  0.0661,  0,0669, 
0.0703,  0.0709,  0.0714,  0.0717, 
0.0707,  0.0699,  0.0688,  0.0675, 
0.0570,  0.0540,  0.0509,  0.0478, 
0.0309,  0.0274,  0.0238,  0.0202, 
0.0033,  0.0004,  0.0001,  0.0001, 
0.0021,  0.0033,  0.0047,  0.0063, 
0.0130,  0.0144,  0>0I57,  0.0167, 
0.0468,  0.0543,  0.0613,  0.0680, 
0.0974,  0.1025,  0.1074,  0.1120, 
0.1299,  0.1323,  0.1345,  0.1363^ 
0.1411,  0.1413,  0.1413,  0.1411, 
0.1359,  0.1340,  0.1320,  0.1297, 
0.1157,  0.1124,  0.1090,  0.1054, 
0.0864,  0.0824,  0.0784,  0.0744, 
0.0539,  0.0498,  0.0456,  0.0415, 
0.0259,  0.0244,  0.0219,  0.0176, 
&END 


0.7090,  0.7238,  0.7385,  0.7533, 
0.8271,  0.8419,  0.8567,  0.8714, 
0.9453,  0.9601,  0.9749,  0.9802, 
1. 0000,  0.9996, 

0.0008,  0.0018,  0.0028,  0,0038, 
0.0256,  0.0291,  0.0323,  0.0353, 
0.0474,  0.0494,  0.0513,  0.0530, 
0.0601,  0.0612,  0.0623,  0.0633, 
0.0676,  0.0683,  0.0690,  0.0697, 
0.0719,  0.0719,  0.0717,  0.0713, 
0.0659,  0.0641,  0.0621,  0.0597, 
0.0446,  0.0413,  0.0379,  0.0344, 
0.0165,  0.0127,  0.0088,  0.0061, 
0.0005,  0.001 1, 

0.0080,  0.0097,  0.0114, 

0.0168,  0.0244,  0.0319,  0.0391, 
0.0743,  0.0804,  0.0863,  0.0920, 
0.1163,  0.1204,  0.1240,  0.1272, 
0.1379,  0.1391,  0.1401,  0.1407, 
0.1406,  0.1399,  0.1389,  0.1375, 
0.1272,  0.1246,  0.1218,  0.1188, 
0.1018,  0.0981,  0.0943,  0.0904, 
0.0704,  0.0663,  0.0622,  0.0580, 
0.0373,  0.0331,  0.0289,  0.0274, 
0.0124,  0.0093, 
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APPENDIX  F:  RVCQ3D  CODE  INPUT 


A  sample  namelist  input  file  for  RVCQ3D  is  shown  below.  Namelist  6 
contains  the  definition  for  a  linear  6%  streamtube  contraction  which  accounts  for  AVDR 
and  the  “Quasi-3D”  effect  in  the  code.  All  solutions  were  run  to  6,000  iterations  and 
consumed,  on  the  average,  only  55  minutes  on  the  NPS  CRAY  (“Sirius”)  and  4  minutes 
on  the  NAS  CRAY  (“Von  Neumann”). 


♦SANGER  CONTROLLED  DIFFUSION  BLADE’ 

4tNLl  M-240*N*49>M[TL=60JMII^^^ 

&NL2  NSTG=44VBT*=J4RS=OTI=^JE»^^  AV2=«., 

AV4=L&END 

^aa*31BClN=14BCEX=ITrMAX=60004RESTI=04RESTO=l,IRES=14CRNT='10000, 

IXRM^dbEND 

4iNL4AMLE=«,24ALLE=43.O»BETE=O4*RAT=O,975O4’0IN=LO,TOIN=LO;G=L4O, 

&END 

ANL5ILT=2^DGE=30JRENR=1.0e6^RNR=.70JPRTR=0.9,TW=0.0,VISPWR*.666666, 
CMUTM«I4»0  &END 

&NLd  OMEGAN),O^LADE=l,NMN-2&eBd 
-.3000  3.000 
1.0000 1.000 
1.0000.9434 
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APPENDIX  G:  DENSITY  RESIDUALS  FOR  CFD  SOLUTIONS 


2000 


3000  4000  feOOO  0000 

lleratlons 


1000  2000 


3000  4ood  feool 

Iterallons 


APPENDIX  H:  Cp  VERSUS  X/C  FORTRAN  PLOTTING  PROGRAM 


cccccccccccccccccccccccccccccccccccccccccccccccc 
C  Program  To  Read  The  Output  From  Rvcq3d.F  And  Grape.F 
C  And  Generate  A  Data  File  For  A  P/PO  Vs  Chord  Plot 
C  Around  The  AIrfolL  Aug.22, 1991 

cccccccccccccccccccccccccccccccccccccccccccccccc 

C  IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  Q(300, 300,4), X(300,300),Y(300,300) 

DIMENSION  U(300,300),V(300,300),P(300,300) 

C 

C  CALCULA  TE  THE  CRITICAL  VELOCITY 
cepe=1.00 
rgas=1.0 
pi=1.0 
ti=1.0 

PRAT=0.685 

rhoi=pi/(rgas*ti) 

ceve=cepe-rgas 

g=1.40 

gP=g+l 

cstar=sqrt(2*g*pi/(gp*rhoi)) 

C 

ISTART=60 

IFINIT=I8I 

C 

READ(7,*)NI,NJ 
print  *,ni,nj 

READ(7,*)((Xa,J),I=l.NI),J=I,NJ), 

&  ((Ya,J),I=l,NI),J=l,NJ) 

READ(7,*)MTL,MIL 

C 

READ(3,*)NI,NJ 

READ(3,*)FSMACH,ALF,RE,TIME 

READ(3,*)(((Qa,J,K),I=I,NI),J=l,NJ),K=I,4) 

GAMMA=I.4 
DO  1234  1=1, MI 
DO  4321  J=1,NJ 
U(I,J)=Q(I,J,2)/QaJ,l) 
va,j)=Qa,J,3)/Qa,j,i) 

Pa,J)=(GAMMA-l)*(QaJ,4)-.5*Qa,J,l)*(U(I,J)**2+Va,J)**2)) 

&  ♦RHOI*CSTAR**2 
4321  CONTINUE 
1234  CONTINUE 
C 

C  DETERMINE  THEXMINAND  XMAX  GRID  POINT  POSITION 

XMIN=0.0 
XMAX=0.0 

DO  2  I=ISTART,IFINIT 
IF(Xa,  1).LE.XMIN)THEN 
IMIN=I 
PRINT  *,IMIN 
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r>  o  o 


XMIN=Xa,l) 

END  IF 

IFPC(1, 1).GE.XMAX)THEN 
IMAX=I 

c  PRINT  *,IMAX 
XMAX=Xa,l) 

END  IF 
2  CONTINUE 
C 

C  FREE  STREAM  STATIC  PRESSURE  USES  THE  XMIN  GRID  POINT 

P0=P(IMIN,NJ)+.5*Q(1MIN,NJ,1)*(U(IMIN,NJ)**2+V(IMIN,NJ)**2) 
&  *RHOI*CSTAR**2 

C 

i=miin 

j=nj 

t0=l/rgas*((g-l)*q(ij,4)/q(ij,l)-(g-l)**2/(2*g)* 

&  (q(ij,2)**2+q(ij,3)**2)/q(ij,l)**2) 
t0=t0*cstar**2 

c  print  *,pO,tO,p(imin,nj),q(imin,nj,l),u(iniin,nj),v(unin,nj) 

C 

C  CALCULATE  THE  MASS  FLOW  RATE  FROM  THE  INPUT  DATA 

C  NORMAUZED  WITH  THE  INLET  AREA 

dniass=pi/(rgas*ti)*(l/prat)**(-l/g)*sqrt(g*(g-I)*rgas*ti/2* 

&  (I-l/(l/prat)**((g-l)/g))) 

c  print  *,’031.  inlet  mass  flow  rate  =',dmass,'*  area' 

C 

(lniasse=q(imin,nj,  l)*RHOI*u(iniin,nj) 
c  print  *,'comp.  inlet  mass  flow  rate  /  area  =',dmasse 
C 

C  PRINT  OUT  DOWN  STREAM  CONDITIONS 

pl=l)(l,l)+.5*q(l,l,l)*(u(l,l)**2+v(l,l)**2)*RHOI*CSTAR**2 

tl=l/rgas*((g-l)*q(l,l,4)/q(l,l,lHg-l)**2/(2*g)* 

&  (q(l,l,2)**2+q(l,l,3)**2)/q(l,l,l)**2) 
tl=tl*CSTAR**2 
c  print  *, 'downstream  condition* 

c  print  *,pl,tl,p(l,l),q(l,l,l),u(l,l),v(i,l) 

C 

C  CALCULATE  THE  MASS  FLOW  RATE  AT  EXIT 

smass=0.0 
ra=0.0 
ba=0.0 

do21  j=l,nj-l 

smass=smass+(q(l  j,  l)+q(l  j+1,  l))*(u(l  J)+u(l,J+l))* 

&  RHOI*CSTAR*(abs(y(lJ+l)-y(lJ)))*0.25 
ra=ra+(q(lj,l)+q(lj+l,l))*(abs(y(lj+l)-y(lj)))*0.5* 

&  RHOI 

ba=ba+abs(y(lj+l)-y(lJ)) 

21  continue 
do  22  j=l,nj-l 

smass=smass+(q(nij,l)+q(niJ+l,l))*(u(niJ)+u(ni,J+l))* 

&  RHOPCSTAR*(abs(y(nij+l)-y(nij)))*0.25 
ra=ra+(q(nlj,l)+q(nij+l,l))*(abs(y(mj+l)-y(nij)))*0.5 
&  *RHOI 

ba=ba+abs(y(ni  j+ 1  )-y(ni  J)) 
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22  continue 
C 

C  AVERAGE  VELOCITY 
va=sniass/ra 
dp=Ta/ba 

c  print  *>01.  downsteani,vd  -  ,va 
c  print  ♦,'ave.  density,dr  =',dr 
C 

C  CALCULATE  THE  INPUT  INLET  VELOCITY 
C  NORMAUZED  BY  CSTAR 
cO=sqrt(g*rgas*ti) 

vin=sqrt(g*(g-l)*rgas*ti/2*(l-(l/prat)**((l-g)/g)))/c0 

c  print  Vinput  inlet  velocity /cO-,vin 

C 

C  CALCULATE  THE  CORRECTED  DOWNSTREAM  VELOCITY 
do23j=nj,l,-l 
vac=u(nij)*vinA^a 
alpha=atan(abs(v(mj)/u(mj))) 
vtc=vac*tan(alpha) 
vc=sqrt(vac**2+vtc**2) 
yy=abs(y(ni,nj)-y(l,nj)) 
c  write(24,*)abs(y(nij)-y(ni,nj))^,vc 

23  continue 
C 

do25  j=2,nj 
vac=u(lJ)*vinA'a 
alpha=atan(abs(v(lj)/u(lj))) 
vtc=vac*tan(alpha) 
vc=sqrt(vac**2+vtc''‘*2) 
c  write(24,*)abs&'(lj)-y(ni,nj))^,vc 
25  continue 
C 
C 

CHORD=ABS(X(IMAX,l)-X(IMIN,l)) 

C 

DO  1  I=ISTART,IFINrr 
DIST=ABS(Xa,l)-X(IMIN,l)) 

XS=DIST/CHORD 

CP=1.0+(-P0+Pai))/(-5*Q(IMIN,NJ,l)*(U(lMIN,NJ)**2 
&  +V(IMIN,NJ)**2)*RHOI*CSTAR**2) 

wcr=sqrt(u(i,  l)**2+v(i,  1)**2) 

C  THE  SURFACE  POINT  Pa,l)  always  EQUAL  0? 

PT=Pai)/P0 

c  WRITE(63,*)XS,WCR 
c  WRITE(64.*)XS,PT 
WRITE(65,*)XS,CP 
1  CONTINUE 
C 

STOP 

END 
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APPENDIX  I:  UPSTREAM  FLOW  ANGLE  FORTRAN  PROGRAM  “angle.r 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  ’an^e.f  READS  FILES  'fort?'  &  'fortS'  WRITES  OUTPUT  TO  FILE  ’fortlO’ 
c  PROGRAM  READS  THE  OUTPUT  FROM  RVCQSDT  AND  GRAPE.f  AND  GENERATES  A 
c  DATA  FILE  CONTAINING  THE  Y-POSITION  AND  INLET  FLOW  ANGLE  AT  A  DESIRED 
c  LOCATION  UPSTREAM  OF  THE  BLADE  IN  TERMS  OF  PERCENT  CHORD.  DATA  IS 
c  CALCULATED  USING  SIMPLE  LINEAR  INTERPOLATION, 
c  SEPTEMBER  29, 1994  WRITTEN  BY  LT.  DREW  WILLIAMS 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


real  beta(300) 

dimension  x(300,300),y(300,300),q(300, 300,4) 
dimension  xx(300),yy(300),u(300,300),v(300,300) 
dimension  uu(300),w(300) 

open(10) 
rewind(10) 
pi=3. 14 15926 

c  Read  in  the  C-grid  created  by  GRAPE 
read  (7,*)  ni,nj 

read  (7,*)  ((x(ij),i=l,ni)J=l,nj), 

&  ((y(ij),i=l,ni)j=l,nj) 

read  (7,*)  mtl,mil 
c  write(10,*)x(mtl,l) 

c  Read  in  the  Solution  file  created  by  R  VCQ3D 
read  (3,*)  ni,nj 

read  (3,*)  fsmach,alplia,reyn,iters 
read  (3,*)  (((q(ij,k),i=l,ni)j=l,nj),k=l,4) 

c  Calculate  the  velocities  u  and  v  at  every  grid  point 
do  100  i  =  l,ni 
do  200  j  =  l,nj 
u(ij)=q(ij,2)/q(ij,l) 
v(U)=q(ij,3)/q(iJ,l) 

200  continue 
100  continue 

c  Define  the  Chord  length  and  then  prompt  for  an  upstream  location 
choi:d=x(mtl,l) 

write(*,*)  Input  Upstream  Chord  Location  in  Percent  (i.e.,  -10,  -20,...)' 
c  (NOTE:  You  must  be  aware  of  the  maximum  upstream  position  available  Jrom  your 
c  GRAPE  output.  Beyond  the  maximum  upstream  location,  the  values  calculated 
c  are  all  the  same  as  that  at  the  maximum.) 
read(*,*)  ax 
axperc=ax/100 
vertx=axperc*chord 
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c 

c 


write(10,*)  chord 
write(10,*)  vertx 


c  Bracket  the  vertical  x-position 

c  (Note:  Depending  on  your  grid  size,  the  7'  index  may  have  to  be  adjusted  to 
c  ensure  the  entire  inlet  flow  region  is  covered.) 
do  300  i  =  1,240 
do400j  =  l,nj 
flag=0 

if  (x(i  j).gt.  vertx)  jsxl=j 
if  (x(ij).It.  vertx)  then 
jsx2=j 
flag=2 
goto  500 
endif 
400  continue 
500  continue 

if  (£lag.gt.l)  then 

c  Calculate  the  u  and  v  velocities  and  compute  the  inlet  angle 
c  Linear  Interpolation  in  the  form  of:  ynew  =yold  +  (dy/dx)*dx 
yy(i)=y(Usxl)+(vertx-x(iJsxl))*0'(iJsxl)-y(iJsx2))/ 

&  (x(ijsxl)-x(ijsx2)) 

uu(i)=q(ijsxl,2)+(vertx-x(ijsxl))*(q(ijsxl,2)-q(ijsx2,2))/ 

&  (x(ijsxl)-x(ijsx2)) 

vv(i)=q(ijsxl,3)+(vertx-x(ijsxl))*(q(ijsxl,3)-q(ijsx2,3))/ 

&  (x(ijsxl)-x(ijsx2)) 

beta(i)=atan(vv(i)/uu(i))*  180/pi 
write(10,*)  5y(i),beta(i) 
endif 

300  continue 
stop 
end 
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